Vol. 74, Part IV, 2004 


ISSN 0369-821 1 




The National Academy of Sciences, India, Allahabad 
%fFr 3T<il>R^, '•TRc!, j^cH^WIcJ 



The National Academy of Sciences, India 

(Registered under Act XXi of 1860) 

Founded 1930 
COUNCIL FOR 2004 
President 

1 Prof Jai Pal Mittal. Ph D (Notre Dame), FNA, FA Sc, FNASc. FTWAS, Mumbai 

Two Past Presidents (including the Immediate Past President) 

2 Prof S K Joshi, D Phi! ,D Sc (he). FNA. FA Sc. FNA Sc. FTWAS. New Delhi 

3 Dr V P Sharma. D Phi! .DSc.FAMS.FESI.FISCD.FNA.FASc.FNASc.FRAS. New 
Delhi 

Vice-Presidents 

4 Dr P K Seth, PhD. FNA. FNA Sc. Lucknow 

5 Prof M Vsjayan, PhD, FNA. FASc. FNA Sc. FTWAS. Bangalore 

Treasurer 

6 Prof S L Srivastava. D Pnil ,F i E T E . F N A Sc , Allahabad 

Foreign Secretary 

7 Dr S E Hasnain, PhD, FNA. FA Sc. FNA Sc. FTWAS. Hyderabad 

General Secretaries 

8 Dr V P Kamboj, Ph D . D Sc . F N A . F N A Sc . Allahabad 

9 Prof Pramod Tandon, Ph D . F N A Sc . Shillong 

Members 

10 Dr Samir Bhattacharya, Ph D . FNA, F A Sc . F N A Sc . Koikata 

11 Prof Suresh Chandra, D Phil , Grad Bnt I R E , FN A Sc . Varanasi 

12 Prof Vtrander Singh Chauhan, Ph D . D Phi! (Oxford). FN A . FN A Sc . New Delhi 

13 Prof Asis Datta. PhD.OSc.FNA.FASc.FNASc.FTWAS .New Delhi 

14 Prof Kasturi Datta. Ph D . FN A . FA Sc . FN A Sc . FT W A S . New Delhi 

15 Prof Sushanta Dattagupta, PhD.FNA.FASc.FNASc.FTWAS. Koikata 
16. Dr Amit Ghosh, Ph D , FNA, F A Sc . F N A Sc , Chandigarh 

17 Prof H S Mam, Ph D (Columbia), FA Sc . FN A Sc , Chennai 

18 Prof G K Mehta. PhD. FNASc. Allahabad 

19 Dr G C Mishra, Ph D . FN A Sc , Pune 

20 Dr Ashok Misra. M S (Chem Engg ), M S (Polymer Sc ), Ph D . FN A Sc .Mumbai 

21 Prof Kambadur Muralidhar, PhD. FNA, FA Sc. FNA Sc .Delhi 

22 Dr Vijayalakshmi Rav.ndranath, Ph D .F N A Sc , FTWAS. Manesar 

23 Prof Ajay Kumar Sood. PhD. FNA. FA Sc. FNA Sc. FTWAS. Bangalore 

Special Invitees 

1 Prof MG K Menon, PhD (Bristol). D Sc (he). FNA. FASc. Hon FNASc FTWAS FRS 
Mem. Pontifical Acad Sc , New Delhi 

2. Dr (Mrs ) Manju Sharma. PhD.FNAAS.FAMi.FlSAB.FNASc.FTWAS. New Delhi 
3 Prof PN Tandon, M S , DSc(h c ), FR C S . FA M S . FN A . FA Sc . FN A Sc . FT W A S . Delhi 
4. Prof Girjesh Govil. PhD. FNA, FA Sc, FTWAS. Mumbai 

The Proceedings of the National Academy of Sciences, India, is published in two Sections Section A 
(Physical Sciences) and Section B (Biological Sciences) Four parts of each section are published annually 
(since 1960) 

The Editorial Board in its work of examining papers received for publication is assisted, in an honorary 
capacity by a large number of distinguished scientists The Academy assumes no responsibility for the statements 
and opinions advanced by the authors The papers must conform strictly to the rules for publication of papers in 
the Proceedings A total 25 reprints is supplied free of cost to the author or authors The authors may ask for a 
reasonable number of additional reprints at cost price, provided they give prior intimation while returning the 
proof. 

Communication regarding contributions for publications in the Proceedings, books for review, subscriptions 
etc. should be sent to the Managing Editor, The National Academy of Sciences. India. 5. Laipatrai Road. 
Allahabad -211002 (India). 

Anr.ja! Subscription for both Sections : Rs. 500.00; for each Section Rs. 250.00; Single Copy : 
Rs. 100.00. Foreign Subscription : (a) for one Section : US $100, (b) for both Sections US $200. 
(Air-Mail charges Included in foreign subscription) 


Co-Sponsored by C S T UP (Lucknow) 



PROCEEDINGS 
OF THE 

NATIONAL ACADEMY OF SCIENCES, INDIA 

2004 

VOL LXXIV SECTION-A PART IV 

A comparative study on the kinetics of oxidation of 
unsaturated acids by quinolinium chlorochromate 

KANCHAN MISHRA, JAI VEER SINGH** and ARCHNA PANDEY* 

Department of Chemistry, Dr. H.S. Gour Vishwavidyalaya, Sagar- 470003, India 
"^Author for correspondence 

"^"^Department of Chemistry, Nehru College, Chhibramu - 209721, Kannauj, India, 

% 

Received October 17, 2001, Revised Apnl 27, 2002, Accepted March 11, 2004 


Abstract 

Kinetics of oxidation of maleic acid and acrylic acid by qmnolimnm chlorochromate (QCC) 
have been investigated m acetic acid - water mixture (50% v/v) at 313 Km presence of 
perchloric acid medium. The order of the reaction is found to be one each with respect to QCC, 
substrate and ¥t for both substrates Increase in the dielectric constant of the medium decreases 
the rate, while vanation in ionic strength has no perceptible change in rate The reaction rates 
have been detenmned at different temperatures and the activation parameters have been 
computed Mechanism consistent with the observed results has been discussed. 

(Keywords oxidation/qxunolmium chlorochromate (QCC) / unsaturated acids / epoxide) 

Introduction 

Many oxidants like Os(VIII)^, Ce(rV)^ and Chloramine-T^ have been used for the 
oxidation of nnsaturated acids. As a part of our earlier mechanistic investigations with 
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unsaturated compounds by QCC'*, we report herein the kinetics of oxidation of 
unsaturated acids by QCC m acetic acid-water (1.1 v/v) solvent in presence of 
perchlonc acid. The mechamsm proposed involves slow attack of protonated QCC on 
the double bond of the substrate in the rate determining step and forms a very unstable 
complex (transient state). These complexes decompose to yield epoxides. 

Materials and Method 

The oxidant, qumolmium chlorochromate^ was synthesised accordmg to reported 
procedure and its purity was checked by iodometric detemunation Solutions of 
maleic acid (Loba) and aciylic acid (B.D.H ) were always freshly prepared. Perchloric 
acid (E Merck) and all other chemicals (A R., B.D H ) were used without further 
purification. 

Kinetic measurements - The reactions were performed under pseudo-first order 
conditions by maintammg a large excess of the substrate over oxidant. The solvent 
used was 50% (v/v) aqueous acetic acid. For both the acids, the progress of the 
reaction was followed by monitoring the decrease m [QCC] at 440 nm usmg 
systronics spectrophotometer. The optical density was measured at various mtervals of 
time. 

Stoichiometry and product analysis - Stoichiometry of the reaction for both the 
acids was found to be IT. Product analysis under kinetic conditions gave only the 
epoxide in each case of study, which was identified by periodate test* 




H 


COOH 


+ Cr02C!0QuH 
COOH 
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COOH O 


+ CrOClOQuH 
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Results and Discussion 

In each case, the oxidation followed a first order kinetics both with respect to 
substrate and oxidant. First order plots of log[QCC] versus time were linear up to at 
least 75% completion of reaction. The pseudo-first order rate constants computed 
from the plots remain unaffected by the change in [QCC] (Table 1 ), establishmg first 
order dependence of the rate on [QCC] in both the cases. At constant [QCC], the rate 
increases steadily with increase in [substrates] (Table 1 ) 

The plots of log kx versus log [substrate] are linear with unit slope (Fig. 1). The 
second order rate constants ki = /:i/[substrate] gave concordant values. The plots of \lk 
agamst l/[substrate] gave a straight line passing through the origin (Fig. 2) This 
confirms first order with respect to substrate and suggests that no mtermediate 
complex is formed between substrate and the oxidant. However, if any complex is 
formed its formation constant would be extremely small. Further UV-visible spectral 
studies did not show any evidence for the formation of QCC-substrate complex (Fig. 3 
&4). 

The direct proportionality between k\ values and [HCIO 4 ] indicate first order 
kinetics with respect to [H^] (Table 1 ). The plots of log kx versus log [H^] gave 
straight line with a slope of unity (Fig. 5). The linear increase in the rate with acidity 
suggests the involvement of protonated Cr(VI) species in the rate determining step, 
similar observations have been made for chromic acid and QFC oxidations’’*. 



Fig. 1- Plot of log kx vs log [unsaturated acid]. [QCC] = 1.88 x 10’’ mol dm \ [HC104]= 1.15 mol dm■^ 
Solvent AcOH- H20= 50% (v/v); Temperature = 313K 
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Fig 2- Plot of [substrate]"*. [QCC] = 1 88x lO'^moldm"^ [HC 104 ]= 1 ISmoldm"^ 
Solvent AcOH H20= 50% (v/v), Temperature = 313K 



Fig 3 - Electronic spectra of QCC in aqueous acetic acid 50% (v/v) without maleic acid (A) QCC m 
aqueous acetic acid 50% (v/v) with maleic acid (B) 
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Added NaC 104 has no appreciable effect on the rate of reaction Similar 
observations were also reported in oxidation of aliphatic alcohols® and 
benzaldehydes*° by QCC. A free radical pathway in the present study is discounted 
This is supported by the absence of polymerization when acrylonitnie was added to 
the reaction mixture. 


Effect of solvent composition - The oxidation of maleic and acrylic acids were 
studied in the solvent containing different proportions of acetic acid and water The 
reaction rate increases with increase in acetic acid content of the reaction mixture. The 
plots of log against HD are linear with positive slope suggestmg an interaction 
between positive ion and a dipole'* This confirms the involvement of protonated 
oxidant species in the rate determining step. 

Effect of temperature - The rate increases with the nse in temperature The 
reaction between unsaturated acids and QCC were studied at four different 
temperatures in the range 303K-318K and the activation parameters were calculated 
(Table 2) from the usual relationships The Arrhenius plots of log k\ versus MT were 
found to be linear (Fig 6). The constancy of the free energies of activation may signal 
the operation of similar mechanism. 


The reaction was foimd to be first order with respect to [QCC], [substrate] and 
[fT]. Moreover the epoxides were the only reaction products observed To explam the 
formation of the product and other observed data, the mechanism in Scheme I is 



WAVEUKNQTH (nmj 


Fig. 4 - Electronic spectra of QCC in aqueous acetic acid 50% (v/v) without acrylic acid (A) QCC ui 
aqueous acetic acid 50% (v/v) with acrylic acid (B). 
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Table 1- Dependence of rate on reactants at 313 K 
Solvent . AcOK - H 2 O 50%(v/v) 


[QCCjxlO’ 

(mol dm'^ ) 

[Maleic acid] x 10^ 

(mol din*^ ) 

[Acrylic acid] x 10^ 

(mol dm'^ ) 

[HCiOJ 

(mol dm'^ ) 

A:iXl0*(s'‘) 

Maleic acid Acrylic acid 

1 13 

- 

2.50 

1 15 

- 

4 32 

1 15 

- 

2 50 

1 15 

- 

4 50 

1 88 

2 50 

2 50 

1 15 

10.5 

4 39 

2 25 

2.50 

2 50 

1.15 

10 7 

4 34 

2 63 

2.50 

2.50 

1 15 

10.6 

4 38 

3 00 

2.50 

2 50 

1 15 

102 

4 41 

3 38 

2.50 

- 

1.15 

104 

- 

3 75 

2.50 

- 

1 15 

10 3 

- 

1 88 

2.75 

2.75 

1 15 

11 4 

4 75 

1 88 

3.00 

3 00 

1.15 

12 4 

5 14 

1 88 

3.25 

3 258 

1 15 

13 3 

5 62 

1 88 

3.50 

3 50 

1 15 

14 2 

6 06 

1 88 

3.75 

3 75 

1 15 

15.3 

6 46 

I 88 

4.00 

- 

1.15 

16.5 

- 

1 88 

2.50 

2.50 

1.38 

12 5 

5 18 

1 88 

2.50 

2 50 

161 

14 9 

5.96 

1 88 

2.50 

2 50 

1 84 

171 

6 77 

1.88 

2 50 

2 50 

2 07 

189 

7.59 

1.88 

2 50 

2.50 

2 30 

21 6 

8.54 
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Table 2- Activation parameters. 


Unsaturated acid E. 


E, 

(KJmor*) 


Scheme I 

AH* 

AG* 

-A 5* 

(KJmol"*) 

(KJmor*) 

(J/Kmor‘) 

52 77 

100.31 

153.12 


Maleic acid 
Acrylic acid 


55.68 


53.10 


102.55 


159 25 
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Fig 5 - Plot of log k vs log [If] [QCC] = 1 88x 10'^ mol dm\ [Unsateated acid] = 2 50x10'^ mol dm‘^ 
Solvent AcOH H2O=50% (v/v). Temperature = 313K 



Fig 6 - Plot of log ki vs 1/T [QCC] - 1 88 x 10'^ mol dm'^ [Unsaturated acid]= 2 50 x 10'^ mol dm‘\ 
[HCIOJ^ 1.15 mol dm'^; Solvent AcOH H20= 50% (v/v). 

The rate determining step may be unstable due to complex formation between the 
protonated QCC and substrate. An electrophilic attack of Cr(VI) being positively 
charged in the protonated QCC on the double bond leads to a four centre intermediate 
state which rearranges to give the ^oxide in die final step. In the proposed 
mechanism (Scheme I) a direct chromium to carbon bond does not account for the 
insensitivity to steric effect very often observed in the oxidation of olefins by 
Cr(VI) Thus the most favourable reaction path may be the tliree centre t5/pe addition 
as depicted in Scheme II. 
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where R= H for acrylic acid and R=COOH for maleic acid. 


Both schemes envisage an oxygen atom transfer from the oxidant. This is in 
accord with the earlier observations made for unsaturated substrates’^'*'*. A closer look 
on the energy of activation of the reactions under investigation reveals that values are 
very close indicating that the electrical and stmctural properties of solvent system are 
little affected by the presence of different substrates. The energy of activation of both 
the acids are nearly equal. This shows that both the acids are oxidised with similar 
rate Consequently the rate law can be proposed as 

[Qcci [H*] (S) 
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Abstract 

Two new highly, sensitive chromogenic spray reagents (AlClrHC104) and (CUCI2/HCIO4) 
have been developed for the detection of some plant growth regulators on silica gel G coated 
chromatoplates The reagents produce different colours m visible light than those in u v light 

(Keywords thm layer chromatography/plant growth regulators/spray reagents) 

Introduction 

Several spray reagents such as 2,4-dinitrophenyl hydrazine’, fenc chloride in 
perchloric acid^ and cenc sulphate^ m concentrated sulphuric acid have been used for 
the detection of indole derivatives on chromatoplates Thomas et alf have used 
solutions of aluminium chloride, thorium chloride, basic lead acetate and neutral lead 
acetate for the detection of flavonoids on paper strips Hramsavljevic et al.^ have used 
cupric chloride for the detection of various oximes. None of the aforesaid reagents is 
ultra sensitive therefore it is worthwhile to develop highly sensitive spray reagents in 
order to analyse the plant growth regulators which are present at a very' low 
concentration in plants Now the efforts are made to develop new and ultrasensitive 
chromogenic spray reagents The results obtamed are described in this paper. 

Materials and Method 

Silica Gel G, perchloric acid, benzene and ethyl acetate (Glaxo, India) and 
aluminium chloride (Danpha, India) were used 

Preparation of chromatoplates 


The plates (15 x 3.5 cm) were prepared by coating a slurry of silica gel G in water 
(2:2 5) with the help of an applicator. The plates were dried at room temperature for 
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4 h, activated at 110°C for 10-15 min , cooled at room temperature and then used for 
analysis 

Preparation of solutions 
Spray reagent I 

Aluinimum chloride (0 5M) was prepared by dissolvmg anhydrous alummium 
chloride (16.63 g) in 250 ml water The spray reagent was prepared by mixing 1 ml of 
this aluminium chloride solution in 50 ml of aqueous perchloric acid (35%). 

Spray reagent 11 

Cupnc chionde (IM) w'as prepared by dissolvmg anhydrous cupric chloride 
(16 82 g) m 125 ml of water. The spray reagent was prepared by mixing 2 ml of this 
cupnc chionde solution with 50 ml of aqueous perchloric acid (35%). 

Test solution 

The solution (1%) of plant growth regulators were prepared in ether or methanol 
The test solution was spotted on TLC plates, the solvent was removed, plates were 
developed m benzene-ethyl acetate as solvent system, dried, sprayed with AICI3- 
HCIO4 and CuCl2-HCI0‘4 separately and then heated m an electnc oven at 120°C for 
1-2 mm to reveal tlie spots. The colour of the spots was noted and then the plates 
were exposed to uv light The colour developed under the uv light w'as also noted 

Results and Discussion 

The Rf values, colour and the lower Imut of detection of plant growth regulators 
by TLC are recorded in Tables 1 and 2 The results show that the spray reagents under 
study are highly sensitive i.e their sensitivities are thousand times more than the 
previously developed ragents®'* (AsCL-AcOH and ASCI 3 -HIO 4 ) for this purpose. 

The advantage of these spray reagents over the previously described spray 
reagents is that they give different colour with different type of plant growth 
regulators, for example mdole derivatives gave light yellow to yellowish violet colour 
with (AICI3-HCIO4) and yellowish brown to reddish brown colour with (CuCL- 
HCIO4) Napthalene acids give blue colour with (AiCl 3 -HC 104 ) and yellow colour 
with (CUCI2-HCIO4) Napthols give pink colour with (AICI3-HCIO4) and grey to 
greyish-purple colour with (CUCI2-HCIO4). Cinnamic acid derivatives like coumaric 
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and protocatechuic acid give brown to green colour with (AICI3-HCIO4) and green 
colour with (CUCI2-HCIO4). However, the spray reagent (AICI3-HCIO4) is more 
sensitive than (CUCI2-HCIO4) as it can detect low levels of plant growth regulators 
(Tables 1 & 2 ). Also the colours developed with spray reagent (AICI3-HCIO4) are 
stable for a longer period than that developed with (CUCI2-HCIO4). 


Table 1- R/ values and colour formation of various plant growth regulators with AICI3-HCIO4 


Compound 

i?/ values m solvent system 

Visible light 

UV light 

(1) 

7.1 

(2) 

6.1 

(3) 

4 1 

(4) 

2 1 

(5) 

Colour 

(6) 

Sensitivity 

(Pg)(7) 

Colour 

(8) 

Sensitivity 

(lig)(9) 

Indole acetic acid 

Oil 

0 12 

0 13 

0.13 

Light 

yellow 

0 25 

Pink 

0 19 

Indole b^tync acid 

021 

0 27 

0.31 

0.35 

Yellowish 

violet 

0 20 

Pink 

0 12 

a-Napthalene acetic 
acid 

0 35 

0 36 

0.39 

041 

Blue 

0 05 

Yellow 

0 12 

P-Napthaiene acetic 
acid 

0 34 

0 36 

0.36 

0 37 

Blue 

0 05 

Black 

0 29 

Phenyl acetic acid 

031 

0 46 

0 54 

0 59 

Orange 

0 15 

Light 

yellow 

0 08 

Phenoxy acetic acid 

0.28 

0 37 

0.40 

0 42 

Pink 

0 30 

Brown 

018 

a-Napthol 

0 76 

0.78 

0 80 

0 80 

Pink 

0 32 

Pink 

0 26 

Kmetin 

0.02 

0 46 

0.55 

0.62 

Black 

0 52 

Brown 

0 39 

Coumanc acid 

016 

0.18 

0.18 

0 20 

Brown 

0 47 

Brown 

0 42 

Protocatechuic acid 

0 07 

0 07 

0.09 

0 09 

Green 

0 35 

Yellow 

0 27 

Vamllic acid 

0.19 

0 20 

0.22 

0.23 

Grey 

0 26 

Black 

0 19 


Solvent system . Benzene . Ethylacetate • Spray reagent ■ AICI 3 -HCIO 4 


The mobile phase such as diisopropylether-dimethylformamide (4:1), n- 
butanolacetic acid-water (5:1:2) and acetone-chloroform-acetic acid-water (8:8.4 1) 
can be used for the separation of plant growth regulators. The best separations are 
achieved in benzene ethyl acetate. 



382 


J C KOHL! et a! 


Table 2~ i?/ values and colour formaiion of vanous plant grov^tfi regulators witli CuCi 2 - HCi 04 


s 

Nc 

Compound 

values m solvent 

s\steni 

Visible light 

U V light 



10 1 

81 

6.1 

Colour 

Sensitivity 

(jig/spot) 

Colour 

Sensitivity 

(pg/spot) 


2 

3 

4 


5 

6 

7 

8 

! 

Indole acetic acid 

0 28 

0 34 

0 41 

Yellowish 

brown 

10 

Light 

green 

05 

2 

Indole buUnc acid 

0 16 

0 24 

0 32 

Reddish 

bro\\ii 

1 5 

Brown 

1 0 

3 

a-Nepthalene 
acetic acid 

0 22 

0 29 

0 39 

Yellow 

1 5 

Pinkish 

brown 

1 5 

4 

p-Napthaiene 
acetic acid 

0 24 

0 32 

0 44 

Yellow 

20 

Piniash 

brom 

1 0 

5 

Kinetin 

0 65 

0 69 

0 74 

Black 

05 

Brown 

05 

6 

Gibberellic acid 

0 69 

0 72 

0.79 

Brown 

1 5 

Blue 

1 0 

7 

Plien>l acetic acid 

0 19 

0 23 

0 29 

Grey 

20 

Blue 

1 0 

8 

Phenoxy acetic 
acid 

0 27 

0 35 

0 46 

Greyish 

purple 

25 

Black 

1 5 

9 

a-Napthol 

0 72 

0 78 

0 82 

Light 

pink 

1 0 

Creaimsh 

yellow 

05 

10 

P-Napthol 

0 74 

0.82 

0 84 

Pinkish 

grey 

20 

Creaimsh 

yellow 

1 0 

11 

Coumanc acid 

0 20 

0 23 

0 27 

Light 

green 

05 

Bluish 

green 

0.25 

12 

Protocatechuic 

acid 

OiO 

0 14 

0 19 

Green 

1.5 

Yellow 

1 0 


Solvent system Benzene . Ethyl acetate (10 ' 1, 8 1 and 6 . 1) 
Spray reagent . CuCla - HCIO4 
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However^ decomposition of indole derivatives was observed in mobile phase 
containing ammonia or formic acid 
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Abstract 

Synthesis of a senes of functional terpolymers, consisting of two electron donating monomers 
namely a-terpmeol and styrene with one electron acceptmg monomer namely acrylonitnle 
(AN), using benzoylperoxide as radical initiator m xylene at 75°C for 2 h have been earned out. 
The system follows ideal kmetics with an overall activation energy 37kJ/mol. The formation of 
the teipolymers is confirmed by the presence of peaks between 7 2-7.6 5 due to phenyl group 
of styrene and at 7.8 - 8 0 8 due to alcoholic group of a-terpineol. NMR spectrum of 
terpolymer has peaks at 120 6 of -CN, 132 5 of -C^Is and 64-66 8 of C-OH. Bands at 3064 
cm“', 2242 cm"' and 3445 cm"* m the FTER spectrum of the terpolymer have mdicated the 
presence of phenyl, cyamde and hydroxy groups respectively The reactivity ratios, determined 
by the Kelen-Ttldos method, are 1 45 for o (a-terpineol + Sty) and 1 3 for rj (AN) The glass 
transition temperature (Jg), determined by DSC, is 59.6“C 

(Keywords ; FTffi/’H/'^C NMR spectra/ reactivity ratios / a-terpineol /styrene / acrylonitrile/ 
terpolymerization / kinetics! 


Introduction 

Despite the large amount of data accumulated for copolymerization of vinyl 
monomers, relatively little synthetic and kinetic information have been published on 
terpolymenzation.*’^ Terpolymerization has become increasingly important from 
commercial view point, because it provides a convenient method for modifying 
polymer properties. 


Styrene (Sty) and acrylonitrile (AN) are highly reactive monomers and can be 
homo, co-and terpolymerized with numerous other vmyl monomers^’'*. These 
monomers frequently copolymerize with isoprene®’’, a-and p-pinene*'’°. However, 
reports on the polymerization of oxygenated derivatives of terpenes except 
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citronellol"''*, geramol’**, linalool’^'® and a-terpineol’''"’* are still scarce A probable 
reason seems to be the inability of the terpenoids to undergo homo-polymenzation 
because of stenc hinderance’®'^ low stabilization energy between monomers and free 
radical transition state^', excessive chain transfer^^ or termination of cyclization^® 
Terpenoid monomers have great relevance and can be used as functional 
poKmers"’’*. Therefore, we have extended our studies to the synthesis and 
characteiization of terpoljmers using monocyclic monoterpenoid, a-terpmeol, styrene 
with acrv'lomtrile We have presented a discussion of the kinetics, thermal properties 
and charactenzation of terpolymers. 


!xoh 

a-terpineol 


a-terpineol, CioHigO, the unsaturated tertiary alcohol, is susceptible to 
polymerization as it contains one double bond 

Materials and Method 

Acrylonitrile, styrene (Merck-Schuchardt) and solvents were purified by the usual 
methods^’^^ a-Terpineol“ (Fluka) was punfied by fractional distillation. 
Benzoylperoxide (BPO) was recrystalhsed twice from chloroform. 

Terpolymenzation 

The terpolymerization kinetics was followed by the dilatometnc technique The 
solution was prepared by taking required quantities of all the three monomers 
alongwith BPO in xylene as an inert solvent. Terpolymerization was carried out at 75 
± 0 rC for two hours under a nitrogen blanket. The terpolymer(s) which precipitated 
from methanol were redissolved m xylene, reprecipitated and dried under vacuum 
The terpolymer formed was treated with acetonitrile to remove polyacrylonitrile, 
cyclohexane to remove polystyrene, and toluene to remove copolymer(s). The rates of 
polymerization (Rp) were calculated from the slope of fte graph between the 
conversion percentage and time. The monomer reactivity ratios were determined by 
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Kelen-Tudos method^^ The FTIR spectra were recorded on a Perkm-Elmer 599B 
spectrophotometer 'H and NMR spectra were recorded on a Varian lOOHA/ 
JEOL 400 LA Spectrometer using CDCl^ as solvent and TMS as internal reference 
Elemental analysis was performed on Perkin-Elmer elemental analyser at 24°C 
Thermogravimetnc analysis (TGA) was performed on a Stanton Redcroft mstrument 
with a heatmg rate of 10“C/mm in nitrogen. Differential Scannmg Calorimetry (DSC) 
was earned out using the DuPont V4 1C Model 2000 analyzer, at a heatmg rate of 
10°C/mm, with a weighmg approximately 10 mg 

Results and Discussion 

The kmetics of terpolymenzation are presented in Tables (1-2) All the reactions 
were associated with a short induction period between 2 and 8 min. The imtiator and 
the monomer(s) exponent values were found to be 0 5 ± 0.01 and unity respectively 

Polymenzation runs were also earned out at 70°C and 85°C to eyaluate the 
activation energy. The Rp is a direct function of temperature and the apparent 
activation energy, determmed from the slope of the Arrhenius plot, its value is 
37kJ/mol (Fig. 1) 

Table 1- Effect of initiator concentration on the rate of terpolymenzation 


Sample 

[BPO] X 10^ 
(moir‘) 

% Conversion 
terpolymer 

RpXlO® 

(moir's-') 

[n] X 10^ 

(dl/g) 

1 

1 83 

62 

3.71 

42 

2 

5 50 

83 

5 88 

40 

3 

9.18 

11 8 

7.58 

37 

4 

12 85 

13 6 

9.54 

3 5 

5 

16 52 

15 5 

9.92 

33 


[AN] 

[Sty] 

[a-trpineol] 
terpolymenzation time 
terpolymenzation temperature 


= 2.34 moir' 
= 1 34 moll-' 
= 0 90 moir’ 
= 2h 

= TSlOffC 
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Table 2- Effect of inonomer(s) concentration on the rate of teipolymerization of a-terpmeol, styrene and 


actylonitnie 


Sample 

a[AN] 

(moir^) 

[Styrene] 

(moll-‘) 

[a-terpmeoi] 

(molH) 

Conversion 

(%) 

RpclO’ 

(moir's"') 

[tllxlO^ 

(dl/g) 

6 

1 00 

1 34 

0 90 

6.80 

5 24 

3 1 

7 

167 

1.34 

0 90 

8.84 

6 45 

35 

3 

2 34 

1 34 

0 90 

11 8 

7 58 

37 

8 

3 01 

1 34 

0 90 

15 98 

8 70 

38 

9 

3 68 

1 34 

0.90 

17 75 

9.77 

4 1 

10 

2 34 

0.57 

0.90 

7.05 

6.02 

25 

11 

2 34 

0 96 

0 90 

9.72 

6 91 

2.8 

12 

2 34 

1 73 

0 90 

17.95 

8 12 

3.4 

14 

2.34 

211 

0 90 

18 95 

8 70 

37 

15 

2 34 

1 34 

0 38 

8 33 

4 16 

29 

16 

2 34 

134 

064 

10 36 

5 01 

2 1 

17 

2.34 

134 

1 16 

14 2 

6 91 

23 

19 

2 34 

1 34 

1 42 

18 1 

7 07 

25 


BPO = 9 18 X 10-^ moll"’ 
teipolymenzation tune = 2 h 


terpolymenzation temperature = 75 ± 0 1°C 

The intrinsic viscosity (T]) of the terpolymers decreases with increase in the 
initiator concentration and increases with increase in the comonomer concentration. 
These observation tend to support the view that although a-terpineol has a hydroxyl 
proton, it might behave also as a chain transfer agent. 

Characterization of terpolymer(s) : 

Solubility : The terpolymer was soluble in aprotic solvents like dimethylforma- 
mide, dioxane, carbon tetrachloride, and chloroform. 
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Fig 1-Arrheiuus plot of rate of polymenzation versus polymenzation temperature, 

[BPO] = 9 18 X 10-^ moir’, [a-teipmeolj = 0 909 moir\ [AN] = 2 34 moir’, 
[Styrene]=l 34 molP , terpolymerization time = 2h 

Elemental Analysis : The percentage of nitrogen in terpolymer has been calculated 
from elemental analysis. The values obtained were 4.0, 4.49, 6.89, 5.61, 11 8, 12 92 
percent respectively (samples 3, 6, 7, 14, 16 and 19) 

FTIR spectroscopy : FTIR spectra were obtained using DMF as a solvent. The 
spectra consists of followmg groups of bands (Fig. 2). 

(i) The range of C-H stretching vibrations due to aromatic ring at 3064 - 3072 
cm''. 

(ii) The range of C-H bending vibrations at 1451-1379 cm"'. 

(iii) The -CsN band of acrylonitrile is present at 2242 cm“'. 

(iv) The -OH band of a-terpineol appears at 3445 cm”'. 
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Fig 2- Fllk spectrum of the teipolymer (sample 3) 
Nuclear Magnetic Resonance Spectroscopy ; 


The chemical shifts of protons, attached to elements other than carbon such as 
-OH, -NH and -SH, to an extent, are influenced by the related phenomena of 
mtermolecular exchange and hydrogen bonding. Signals appear in the NMR spectra, 
due to -OH (hydroxyl protons) with species of small molecular weight where 
mtermolecular association is not hindered, in the region 5-3-S.5 ppm (hydroxyl 
proton of methanol appears at 5=3.3 ppm, whereas fliose of ethanol appears at 5=5 4 
ppm)^*’^®. However, with many large molecules, hydroxyl protons often appear near 
5=8 ppm even at relatively high concentration, partially due to steric effects and 
partially due to resonance stabilization. Therefore, the peaks of -OH group appears 
in the range of 7.8 6 to 8.0 5 in the NMR spectra of a-terpineol as well as in 
terpolymer. 
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The NMR spectra of pure a-terpineoI(I) and those of terpolymers (II) (Fig. 3) 
show the following peaks which were assigned (IH, OH) singlet at 7.7-8 0 5, (3H, 
CHj) singlet at 0 9-1 0 5, (2H, CH^) doublet at 4.2-4.7 6, (IH, CH) triplet at 1 8-2 2 5. 
However, an additional peak in terpolymer is assigned for (SHjCeHs) doublet at 7.2- 
765 




Fig. 3 (H) 

Fig. 3 - 'H NMR Spectnim of the terpolymer (sample 3). 


6 (ppm) 
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Fig 4~ Kelen-TMos plot of determination of reactivity ratios 

'^C-NMR 


The spectrum of the terpolymer (Fig. 5) shows the following peaks (-CH^), 
21ppm, (-CHj), 42ppm, (-CH) 52-5 3ppm, (-C^Hj), 132ppm; tertiary C-OH 64- 
66ppm, (-CN), 120ppm. 



Fig. 5- NMR spectrum of teqiol)iner (sample 3). 
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The bands and peaks in the FTIR, ’H NMR and NMR spectra respectively, 
confirm the presence of styrene, acrylonitrile and a-terpmeol units in the terpolymer 
chain 

Differential Scanning Calorimetry : 

The DSC scan of terpolymer (Fig. 6) showed a well pronounced endothermic 
transition around the temperature 59°C. The endothermic transition in the terpolymer 
can be attributed to the formation reaction as reported by various authors for 
polystyrene, polyMMA, and its terpolymers^®’^* with terpenoids’®. The values of glass 
transition temperature onset temperature (T), and the peak temperature (T ) are 
as follows . 

T = 59°C; T = 142.8; T = 145.6, T -T = 2.8 

g ’ O P ^ p O 



D 50 100 150 

Temperature (OC) 


Fig. 6- DSC curve of the terpolymer (sample 3) 
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The onset temperature (To) is obtamed at the intercept of the tangents to the base 
Ime at the lower temperature side of the endothermic peak. T? is the temperature at 
which the bulk of the terpolymer has undergone a dehydration reaction and the 
difference between Tp and To (T^-To) is the measure of the overall rate of reaction The 
smaller the difference the greater is the rate of reaction. 

Thermogrcmmetnc Analysis 


Results from thermogravimetric analysis (Fig. 7) show that the terpolymer is 
stable upto 370°C and rapid decomposition starts at 390°C 



Terpolymer Composition and values oj Reactivity Ratios : 

The moles of [ct— teipineol] and [Sty] formed have been calculated from the peak 
area of hydroxy(7 8—8 0 8) and phenyl protons (7.2— 7.6 8). The moles of acrylomtrile 
formed have been calculated from elemental analysis (%N). The Kelen-Tiidos 
approach has been used for the evaluation of and r^ for the monomer(s) as follows : 

n = M ^ (1-^) 

a 


Deriv. Weight (%/mln) 
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where t| = GI(a-^H) and ^ = Hl{a. + H). The transformed variables G and H are given 
by 

^ _ {M,V[M,]{d[M,Vd[M,]-\) 
d{M,Vd[M^] 

d[M,]ld[M2] 

where, M\ and Mi are moles of styrene + a-terpeniol + AN 

We have calculated the parameter a by taking the square root of the product of the 
lowest and highest values of H for the terpolymenzation series. Graphical values of 
Sty+a-terpineol/ AN yields /•,=1.45 and r,=l 3 (Fig 4) (Table-3) The values of n and 
ri favour random terpolymerization^^’^^. 


Table 3- Composition ofterpolymers 


Sample 

Monomer feed 

[(Sty+a-terp )/AN] 

AN 

Mole fraction 

[Sty] 

[a-terp ] 

n 

4 

6 

2 22 

0 08 

0.58 

0 33 

2 98 

0.65 

13 

128 

0 09 

0 58 

0 34 

2 90 

0.39 

3 

0.95 

0 12 

0 51 

0.36 

2 28 

0.35 

16 

106 

01 

0 48 

0 41 

2 42 

0.35 

8 

0 97 

0.19 

0.41 

0 38 

1 60 

0.48 

9 

0.74 

0.20 

0.42 

0 36 

1.44 

0.36 


Conclusion 


' The radical terpolymerization of [AN] with [Sty+ a-terpineol] results in 
functional terpolymer. 
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Abstract 

A computer assisted investigation has been carried out on the nature of ternary complexes of 
uranium(VI) with oxalic acid as primary ligand and malomc acid as secondary ligand Xhe 
approximate formation constants have been calculated with SCPHD from the pH metric data 
The formation constants thus obtained are refmed with MINIQUAD 75 The selection of the 
best-fit chemical model is based on the statistical parameters and residual analysis The major 
complexes fonned are : [U02(C204XH3C304)2]^', [U02(HC204)(H3C304)] and [U02(C204)2 
(H3C304)]^' The distribution patterns of the different species with the pH values showed that 
[U02(C204)(H3C304)2]^' IS the predominant species The nonlmear trend of stabilities of the 
complexes with surfactant concentration is explamed based on the dielectric constant of the 
medium. Stem layer and the electrostatic mteractions 

(Keywords formation constant/temary complexes/micelles/uranyl lon/oxalate/malonate) 

Introduction 

Ternary complexes of uranyl ion with biologically important ligands^"^ and with 
some carboxylic acids have been stiidied^’^, Uranyl complexes of dicarboxylic acid 
have also been studied in aqueous^"^ and in aqua-dioxane mixtures^^’^\ The stability 
constants of the ternary complexes of copper(II) with malonic acid as the primary 
ligand and oxalic, maleic, fumaric, succinic, phthalic and salicylic acids, 
pyrocatechol(0,0 donors), 0-aminophenol, 8-quinolinol (0,N donors) and 2,2- 
bipyridyl, ethylene diamine(N,N donors) as secondary ligands, have been reported^^. 
Ternary complexes of uranyl ion with 3,5-dmitrosalicylic acid and other carboxylic 
acids were also reported*^. A considerable number of complexes of - 
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carboxylate with aimnoacids*''''’ have been r'^ported but there are only few reports^-''' 
on U02-dicarboxylate with other dicarboxylic acids. Oxalate ion is also used for the 
separation of Pu(IV) and Th(IV) from U(VI) by precipitation**. 

Since fatty acids are useful reagents for the extraction of uranyl ion, a study on the 
stabilities of uranyl complexes with dicarboxylic acids, assumes significance. Hence, 
protonation’^ and uranyl complex eqmlibria of malonic^“ and oxalic acids^* have been 
studied In the present paper the stability of ternary complexes has been reported In 
these studies micellar medium is chosen for two reasons . 1) micelles decrease the 
dielectnc constant of the medium, which is essential in the extraction of uranium and 
2) micelles mimic the biological membrane and compartmentalize the chemical 
reactions. 


Materials and Method 

Solutions of uranyl nitrate, oxalic acid and malomc acid were prepared by 
dissolving G R grade (E. Merck, Germany) samples in triple distilled water Aqueous 
solutions of AR grade sodium laurylsulphate (SLS) (Qualigen), cetyltrimethyla- 
mmonium bromide (CTAB) (BDH) and TritonX-100 (TX) (E. Merck) were ^so 
prepared in triple distilled water. Alkalimetnc titrations were carried out in the 
medium containing varying compositions of SLS (0.5 - 2.0%w/v), CTAB (0.5 - 
2 5%w/v) and TX (0.5 - 2 5%v/v) maintaining an ionic strength of 0.17mol dm‘* with 
sodium chlonde at 303.0±0.1K temperature using constant temperature bath. 
Systromcs (model 335) pH meter was used to monitor changes in hydrogen ion 
concentration. The glass electrode was equilibrated m a well stirred micellar solution 
containing mert electrolyte. The effect of variations in asymmetry, liquid junction 
potential, activity coefficient, sodium ion error and/or dissolved carbon dioxide on the 
response of glass electrode was taken into account in the form of correction factor (log 
F). The log F was computed from the simulated acid-base titration data (PHCi) by 
SCPHD program^^ for each of the solvent compositions. Its value was used to convert 
pH meter dial readmg into logarithm of reciprocal of hydrogen ion concentration 
(pHEi) according to the eqn.(l). 

log F = PHCi - PHEi (1) 

Strong acid was titrated with alkali at regular intervals to check equilibration of 
the glass electrode The calomel electrode was refilled with micellar solution (only TX 
and CTAB but not SLS since it forms precipitate with KCl) of equivalent composition 
as that of the titrand. In each titration, the titrand consists of approximately 1 mmol 
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mineral acid (HCIO4 in case of SLS and TX and HCl in case of CTAB because CT 
AB forms precipitate with HCIO4) in a total volume of 50cm^. Titrations of solutions 
contammg different ratios (1:2.4 and 1:4 2) of metal to oxalic acid to malonic acid 
were earned oi-t with 0.2 mol dm'^ sodium hydroxide. 

Modelling strategy 

A chemical model for complex equilibria consists of stoichiometric coefficients 
and logarithm of stability constants. TTie best-fit chemical models were amved at by 
adopting a modelling strategy Various models were obtained from the experimental 
data usmg MINIQUAD75^''. Some heuristics were followed^^ in the refinement of the 
stability constants The best model among them was chosen based on statistical 
parameters. 


Results and Discussion 


Complex equilibria 

The coordination number of uranyl ion is six^^ and three bidentate ligands are 
sufficient to satisfy the coordination number of the metal ion Hence, the total number 
of primary' and secondary ligands together is restricted to a maximum of three in 
generatmg the possible ternary species for modelling. Several thumb rules based on 
chemical principles and well established practices in multiple linear regression 
analysis^^ were followed in the refinement process which have resulted in the final 
model consisting of the species MLX 2 H , MLXH 2 and ML 2 XH^' The best fit model 
was chosen as that with the lowest standard delation (SD) m formation constants and 
minimum U^on (corrected for degrees of freedom) which was corroborated by other 
statistical parameters^® like ciy'Stallographic R etc., which are given in Table 1. The 
values of skewness between -0.06 and 0.91 evince that the residuals form a part of 
normal distribution. Kurtosis is the degree of peakedness of a distribution, usually 
taken relative to a normal distnbution. If this value is equal to 3, the distribution of 
residuals is normal (mesokurtic) i.e. not very peaked or flat-topped. A high value 
indicates that the distribution is peaked (lepto^rtic) and a low value represents a flat 
topped distribution called platykurtic. 

The effects of errors in the influential parameters'^ like the concentration of 
ingredients on the stability constants were studied by incorporating errors in them in 
the range -1 to + 1 %. The results show that the order of affecting the magnitudes of 
stability constants is due to errors in the concentrations of alkali> acid> ligandZ > 
ligand 1 > metal. The increased SD in stability constants on the introduction of errors 
confirms the correctness of the proposed models. 



402 


BB,V SAILAJAe/^/. 


Table 1- Best fit chemical models of mixed-ligand complexes of Uranmm(VI) with oxalic and malomc 
acids m water-surfactant mixtures 

Temp = 303K p= 0 1 7 mol dm*^ pH range = 2 0-6 0 


log Pxnlxh(SD) 


% surfactant 

MLX 2 H 2 

MLXH 2 

ML 2 XH 

NP 

Hcorr 


R 

Skewness 

Kurtosis 

00 

23 44(8) 

20 04(1) 

18 06(1) 

93 

6 85 

20 37 

0 0114 

-0 06 

3 15 

SLS(%wyv) 










05 

22 17(2) 

19 11(2) 

17 30(5) 

94 

13 17 

26 82 

0 0152 

0 39 

2 81 

10 

23 05(1) 

20 29(1) 

18 16(3) 

88 

9 91 

14 55 

0 0135 

0 42 

2 83 

1 5 

24 33(1) 

20 35(7) 

20 12(2) 

82 

2 28 

2 29 

0 0066 

-0 15 

2 90 

20 

23 66(2) 

20 79(2) 

19 21(3) 

76 

5.96 

29 16 

0 0109 

-0 10 

5 06 

TX(%v/v) 










05 

24 54(8) 

20 50(9) 

20 06(1) 

94 

1 39 

16 63 

0 0059 

-0 15 

3 14 

1 0 

23 40(5) 

21 69(3) 

19 12(6) 

93 

32 37 

51 67 

0 0249 

0 91 

3 44 

I 5 

22 70(1) 

19 73(1) 

18 06(2) 

94 

13.74 

70 21 

0.0171 

0 59 

3 02 

20 

22 29(1) 

19 76(1) 

17 58(2) 

91 

12 16 

70 02 

0.0168 

0 79 

3 42 

25 

23 42(2) 

20 42(2) 

18 50(4) 

95 

11.96 

48 93 

0 0157 

0 57 

2 89 

CTAB(%w/v) 









05 

23 18(1) 

20.33(2) 

18.19(4) 

93 

13.44 

47 31 

0 0172 

0 69 

3 24 

10 

23 71(2) 

20 77(2) 

18 76(5) 

90 

10.57 

37 96 

0 0156 

0 63 

3 13 

15 

24 29(7) 

20 33(8) 

19 15(3) 

90 

1 85 

23 56 

0 0066 

-0 35 

4 19 

20 

22 74(1) 

20 08(2) 

17.92(4) 

95 

27 28 

96 84 

0 0239 

0 86 

3 72 

25 

23.65(2) 

19 80(2) 

18 46(7) 

91 

3 23 

16 92 

0.0085 

-0 43 

4.33 
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Distribution diagrams 

The concentrations of binary metal-ligand complexes were found to be less than 
10% in the presence of ternary complexes Hence they are not presented in the 
distribution diagrams (Fig. 1). Smce ternary complexes are more stable than binary 
complexes, the latter are converted to the ternary complexes in the presence of 
secondary ligand. The change in the stability of ternary complexes as compared to 
their binary analogues was quantified^^’^®. 



1 3 5 7 9 

pH 

Fig 1- Distribution diagrams of various mixed ligand species. FL FX * concentrations 
of protonated primary and secondary ligands, respectively 
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In one of the approaches the difference in stability (AlogAO for the two reactions viz. 
ML with X and M(aq) with L and X (where X is a secondary' ligand) is compared with 
that calculated purely on statistical grounds (eqn 2) 


A log K = log Xjtfec — loE — log K^x ~ log (2) 


Both electrostatic theories of bmary complex formation and statistical arguments 
suggest that more coordination positions of a given multivalent hydrated metal ion 
will be available for the first ligand L.an for the second Hence, Ae usual order of 

stabilitN' bolds good. This leads to a negative value for AlogAT The 
statistical values of AlogAT for bidentate L and X are -0.4, -0 6 and between -0 9 and 
-0 3 for octahedral, square planar and distorted octahedral complexes, respectively 
■Whenever the experimental value of AlogAT exceeds the statistical value, it can be 
inferred that the ternary complex is formed as a result of mteraction of MLHh with 
MXHh Negative values of AlogX can be understood as the secondary ligand forms a 
more stable complex with hydrated metal ion than with ML, which does not mean that 
the teraarj' complex is absent. The ternary complexes with N,N donor pnmary ligands 
are observed^®'^^ to have positive AlogX values if the secondary ligands are O, O. 
donors, negative if the secondary ligands are N,N donors and midway between the 
tw o if the secondary ligands are N,0 donors 


Another approach to quantify the stability of ternary complexes was based on the 
disproportionation constant (logX) given by eqn. 3, which corresponds to the 
equilibrium given by eqn. 4. Under these equilibrium conditions, one would expect 
50% ternary complex and 25% each of the binary complexes to be formed and the 
value^^® of logJ to be 0 6 Any value greater than 0.6 accounts for the extra stability 
of MLX, which results from several factors like charge neutralization, chelate ring 
size, stacking interactions, hydrogen bonding and charge transfer interactions 


log X - 21og AT^ - log K ^2 ~ log X ^2 


(3) 


ML2 + MX2 ^ 


^ 2MLX 


(4) 
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The values of Alog^ and log.Y could not be calculated for the present systems, 
using the transformations given in Table 2 due to the absence of the relevant binary 
complexes 

A perusal of Fig. 1 reveals that at very’ low pH the concentrations of MLX 2 H 2 "' 
and ML 2 XH^' complexes are less than those of protonated ligands. As the pH 
increases the concentrations of MLX 2 H 2 ^’ and MLaXH^' go through maximum but the 
concentration of MLXH 2 species is high at lower pH and gradually decreases as the 
pH increases. In the pH region 1.5-3 .5 oxalic and malonic acids exist as LH (HC2O4) 
and XH 2 (H4C3O4) respectively'^. These protonated ligands interact with the metal ion 
to form MLXH 2 . Since the pKa values of oxalic acid are lower than those of malonic 
acid, oxalic acid is readily deprotonated. Hence the proton in protonated ternary 
complexes must be present on malonate but not on oxalate. The possible formation of 
other ternary species is shown in the following equilibria. 


+ LH2 + XH2 MLXH2 + 2H^ (5) 

(or) 

MLH^ + XH2 MLXH2 + H^ (6) 

MLXH2 -O’ MLXH' + H^ (7) 

MLXH’ + LH“ ML2XH^' + H^ (8) 

(or) 

ML + XH2 MLjXH^' + H^ (9) 

ML + 2XH2 <-> ML2XH^- + 2H^ (10) 


The possible structures of these temaiy complexes are given in Fig. 2. The 
structures of the complexes are proposed based on the literature and chemical 
principles. Literature also shows that the coordination number of uranyl ion is 6 and 
oxalate and malonate act as bidentate ligands. It also shows that water occupies one of 
the coordination positions. 
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The magnitudes of protonation constants of malonic acid are higher than those of 
oxalic acid'^. Therefore, mono- or di- protonated malonate ion exists in the pH range 
where deprotonated oxalate ion can exist. For example, MLXH2 (L = oxalate), 
(X=malonate) can be represented as M(LH2)(X), ML(XH2) or M(LH)(XH) and the 
possibility of M(LH2)(X) can be ruled out because oxalic acid cannot exist in the 
protonated form when malonic acid is present as deprotonated ligand. The possibility 
of ML(XH2) can also be ruled out since diprotonated ligand (XH2) cannot form the 
complex due to unavailability of coordinating atoms. Hence, the MLXH2 species must 
be present as M(LH)(XH). Monoprotonated oxalic or malonic acid shall always act as 
monodentate ligand due to the availability of one coordinating atom. Thus the 
structure of MLXH2 shall be as represented in Fig. 2 . 

Similarly [MLXjHj)]^' can be represented as [ML(XH)2]^’ and [ML2XH]^- as 
[ML2(XH)]^’ whose structures are also s. 3 wn in Fig. 2 . 

Table 2- Calculation of Alog/f and logA" from overall stability constants 


Aiog^ii }2 “ ^t>gPlll 2 * '^ogPioio 

~ ^ogPm 2 * iogpiioi ‘logpioii 

~ iogpin 2 - logPnoo 'logpio !2 

AlOgA:H 22 = l 0 gpn 22 • iogP ,,02 'lOgPi 020 

= 10 gPn 22 * iOgPuOl -Iogpl 021 

~ JogP!l22 • iOgPllOO •'i0gpi022 

AIog/iri 2 i j ” iogpnii - logp| 20 i “iogpioio 

^t)gP| 2 ii - iogp 1200 -Jogp 1011 


^OgA"in2 “2 10gpin2 • ^t>gpi204 ’logpl020 
=2 iOgPiij 2 - logpi203 -iogPi 021 
=2 IogPiil 2 ' ^Ogpi 202 •'^OgPl 022 
-2 logpiii2 - ^OgPi201 '•logPi023 
=2 IogPin2 *■ ^OgPl200 ■■iOgPl024 
logPii22 - Iogpi204 "iOgPl040 

=2 10gPn22 ' iOgPi203 '^OgPl041 

=2 lOgpji22 - iOgPl202 'Jt)gpi042 
=2 l0gp| 322 • ^Ogpl201 "^OgPl043 
=2 10gpii22 “ lOgPj200 ■'iOgPi044 
^OgA"j 2 H =2 logpi 2 n » iogPi40i *IogPi 02 i 
~2 iOgPi2ii - Iogpi4oo -logPi022 
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MLX2H2 



MLXH2 



OCOCI^OOH 

ML2XH 


Fig. 2- Probable structures of uranyl-oxalate-malonate chelates - 



%surfactant 


%surfactant 




% surfactant 


■ Variation of fonnation constant (log with percentage concentration of surfactant. 
(A) CTAB, (B) TX and (C) SLS. ©MLXjHz HMLXKj ▼ MLjXH 
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Effect of micelles 

The vanation of stability constants (logp) with percentage of surfactant is shown 
in Fig 3 Siirilar to binary complexes*^’^°, the stabilities of ternary complexes also 
exhibit nonlinear trend This is also reflected in the maximum concentration of the 
species and the pH at maximum percentage concentration of the species (Table 3). 
The dielectnc constant of the medium decreases^’’^* with increasing concentration of 
the surfactant. The variation in the stability of the charged complex species depends 
on the polarity of the medium, charge on the micellar surface (Stem layer), and on 
electrostatic attraction or repulsive forces operating between the complex species and 
the micellar surface. The species [MLX 2 H 2 P' and [Ml^XH]^' should be stabilised by 
CTAB due to electrostatic interactions but they should be destabilised due to the 
decreased dielectric constant The magnitudes of these opposing factors determine the 
trend in variation of stability constant with increase in CTAB content as shown in 
Fig 3A The stability of MLXH 2 is little affected by the medium because it is a neutral 
species Due to similar reasons the stability of the ternary complexes varies 
nonlinearly with the increased SLS content (Fig. 3C). As TX is a neutral surfactant, 
the stabilities are decreased due to the decreased (Fig 3B) dielectric constant and 
decreased electrostatic interactions. 


Table 3- Effect of surfactant on the pH at vluch maximum concentration of the species is observed for 
the ternary complexes of uramum(VI) with oxalic and malomc acids 


Surfactant 


pH(%concentration of species) 



MLX 2 H 2 

MLXH 2 

ML 2 XH 

00 

3.59(83 1) 

2.0(73.3) 

5 39(26 4) 

CTAB(%w/y) 




05 

4.33(73 1) 

2.0(83 3) 

5.72(25.1) 

1 0 

4 33(75.3) 

2.0(83 8) 

5 79(27.8) 

1.5 

4.45(78 9) 

2.0(84.1) 

5.69(28.9) 

20 

4.43(64 9) 

2 3(82 1) 

5.57(23 9) 

25 

3 47(93 1) 

2.0(67.4) 

5.56(14.1) 


Table 3 Contd... 
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Table 3 Contd 


TX(%v/v) 


05 

3 53(94 1) 

2 0(56 8) 

5 83(46 6) 

1 0 

4 63(92 8) 

2 93(98.2) 

5 95(41.7) 

1 5 

4 35(75 4) 

2 08(83 6) 

5 65(34 1) 

20 

4 34(60 9) 

2.11(86 7) 

5 57(27 2) 

25 

4 35(75 9) 

2 03(83 9) 

5 74(29 5) 

SLS(%w/v) 




05 

4 21 (78 4) 

2.05(81 7) 

5.70(25 5) 

1 0 

4 20(71 9) 

2.10(83 3) 

5 71(22 0) 

1 5 

4 10(95 1) 

2 21(84.2) 

5 33(37 1) 

20 

4.50(80 2) 

2 02(977) 

6 00(26 2) 


Conclusions 

(1) The interaction of oxalic and malonic acids with uranyl ion in the pH range 
2“6 forms 

[U02( C204)(H3C304)2]^ [U02(HC204)(H3C304)] and [U02(C204)2(H3C304f 
species. 

(2) The concentrations of binary complexes of uranium with oxalic or malonic 
acid are very small (less than 10%) compared to those of the ternary 
complexes. 

(3) The magnitudes of stability constants of the ternary complexes exhibit non- 
linear rend with the concentrations of the surfactant This trend is due to the 
decrease in the polarity of the solvent with increased surfactant 
concentration, which destabilizes the charged complex species. The second 
factor IS the electrostatic interactions between the complex species rnd the 
charged Stem layer of the micelles. 
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Abstract 

The present paper denves semi-circle theorems, for magnetorotatory double-diffusive 
convection problems of Veroms and Stem types, that prescribe upper limits for the complex 
growth rate of oscillatory motions of neutral or growing amplitude m such a manner that it 
naturally culrmnates m sufficient conditions precluding the non-existence of such motions for 
an mitially bottom heavy as well as an initially top heavy configurations. Further, results for 
double-diffusive convection problems with or without the individual effects of a rotation or a 
magnetic field follow as a consequence 

(Keyvrords' magnetorotatory/double-diffusive/convection/oscillatory/growth rate) 

Introduction 

Double-diffiisive convection, with its archetypal case of heat and salt, generally 
referred to as thermohaline convection, has been extensively studied in the recent past 
on account of its interesting complexities as w'ell as its direct relevance in many 
problems of practical interest in the fields of the Limnology, Oceanography, 
Geophysics, Astrophysics and Chemical Engineering etc. For a broad and a recent 
view of the subject one may be referred to Brandt and Fernando’. Two fimdamental 
configurations have been studied in the context of thermohaline instability problem, 
the first one by Stem^ wherein the temperature gradient is stabilizing and the 
concentration gradient is destabilizing and the second one by Veronis^ wherein the 
temperature gradient is destabilizing and the concentration gradient is stabilizing. The 
main results derived by Stem and Veronis for their respective configurations are that 
both allow the occurr^ce of a stationar>' pattern of motions or oscillatoiy motions of 
growing amplitude provided the destabilizing concentration gradient or the 
temperature gradient is sufficiently large. However, stationary pattern of motion is the 
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preferred mode of setting in of instability in case of Stems configuration whereas 
oscillatory motions of growmg amplitude rje preferred m Veroms' configuration 
Further, these results are independent of the initially gravitationally stable or unstable 
character of the two configurations. It is important to note that Veronis' work is 
restncted to dynamically free boundaries whereas Stem's work assumes the 'principle 
of exchange of stabilities'. Keeping in view the foregoing discussion, thermohaline 
configurations of Veronis and Stem types can therefore be further classified into the 
following two classes; 

(i) The first class in which thermohaline instability manifests itself when the 
total density field is initially bottom heavy and 


(ii) The second class in which thermohaline instability manifests itself when the 
total density field is initially top heavy. 


Baneqee and co-author|^ have derived a characterization theorem for rotatory 
(hydromagnetic) double-diffusive problem of Veronis' type that disproves the 
existence of oscillatory motion of growing amplitude in an iratially bottom heavy 


configuration if (i) x < 0 (xo] < 0) and (ii) 31 ^ < 


11 %' 

4 



where x, 0, 01 and 


9I5 respectively denote the Lewis number, the thermal Prandtl number, the magnetic 
Prandtl number and the concentration Rayleigh number. The restriction (i) in the 
above results of Baneijee and co-authors may be physically justifiable in certain 
situations, however, it is not mathematically palpable. Further, the sufficient character 
of condition (ii) coupled with the nature of their mathematical analysis strongly 
suggest the possibility of the derivation of an upper bound for the modulus of the 
complex growth rate of an arbitrary oscillatory perturbation which may be neutral or 
unstable that will automatically take care of condition (ii) and yield the results of 
Baneijee and co-authors without the restriction (i) and will also be uniformly 
applicable for initially bottom heavy as well as initially top heavy configurations. 


Motivated by these considerations, the present paper purports to deal with the 
more general magnetorotatory double-diffusive convection problems of Veronis' and 
Stem's types and derives semi-circle theorems that prescribe upper limits for the 
complex growth rate of oscillatory motions of neutral or growing amplitude in such a 
manner that it naturally culminates in sufficient conditions precluding the non- 
existence of such motions for initially bottom heavy as well as initially top heavy 
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configurations. Further, results for simple double-diffusive problems with or without 
the mdividual effects of a rotation or a magnetic field follow as a consequence. 

Basic Equations and Boundary Conditions 

Following the usual steps of linear stability theory the non-dimensional linearized 
perturbation equations governing the magnetorotatory double-diffusive convection 
problem with a uniform rotation and a uniform magnetic field both acting m the 
vertical direction opposite to the force field of gravity are given by (c f Veroms^ and 
Chandrasekhar®). 




cr; 


w = 9{a^Q-9i,a^-QD(D^)h, + TDC, 


2 ^ 2 \ 


- p) 0 = -w 


r 


\ 


<1) = . 


W 


\ cr / 


h, = -Dw 


f 


\ ct; 


C, = -QD^ - Dw 


-a^ = 

\ CT ; 

together with the boundary conditions 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


w = 0 = 9 = <i) = Dw - hi = l^ = D^ at z = 0 and z = 1 
(both boundaries rigid and perfectly conducting) 


( 7 ) 
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In the above equations, z is the real independent variable such that Q<z<\,D = didz is 
the differentiation with respect to z; o, t and Oi are positive constants; g and T are 
non-negative constants; SI and 91, are also constants that could be positive or negative; 
p=Pr + ip, IS a complex constant in general and the dependent vanables w, 6, cj), h,, 
and ^ are complex-valued functions of the real variable z The meanmg of the 
svinbols from the physical point of view are as follows ■ z is t!ie vertical co-ordmate, 
d/dz is the differentiation along the vertical direction; is the square of the wave 
number; a is the thermal Prandtl number; t is the Lewis number, ai is the magnetic 
Prandtl number, 91 is the thermal Rayleigh number; 91, is the concentration Rayleigh 
number, Q is the Chandrasekhar number; T is the Taylor nuinber; p is the complex 
growth rate, w, /j,, ^ and ^ are the vertical components of velocity, magnetic field, 
v'orUcity and current density respectively; 6 is the temperature and <{> is the 
concentration. 


Mathematical Analysis 

System of eqn. (1) - (6) together with the boundary conditions (7) constitute an 
eigenvalue problem for the complex growth rate p = Pr+ tp, for given yalues of the 
other parameters, namely, a^, a, 91, 9I„ Q, T, Oi and t and a given state of the system 
is stable, neutral or unstable according as pr is negative, zero or positive Further, 

(a) p, 7 ^ 0 and pr > 0 describe oscillatory motions of neutral or growing 
amplitude; 

(b) 91 > 0 and 91, > 0 and either Q = 0 = TotQ=^QotT=Q respectively 
describe Veronis' thermohaline configuration (VTC) or rotatory VTC or 
hydromagnetic VTC; 


(c) 91 < 0 and 91, < 0 and either Q = 0 = TotQ = 0otT=0 respectively 
describe Stem's thermohaline configuration (STC) or rotatory STC or 
hydromagnetic STC; 

9? 

(d) X = — <1 

91 , 

and 

(e) X. > 1 (x > 1 jdescribes an initially top heavy VTC (STC). 


91 


describes an initially bottom heavy VTC (STC); 
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Finally, if /?r - 0 => /j, - 0, Va^, then the principle of exchange of stabilities is valid, 
otherwise we have overstability at least when instability sets in as certain modes. 

We prove the following Theorems ; 


Theorem I (A Semi-Circle Theorem for Magnetorqtatory VTC) 


If (p, w, 0, (j), hz, Q ^), p=Pr+ ipi, ^ 0, A ^ 0, is a nontrivial solution of the eqn. 
(1) - (6) together with the boundary conditions (7), and 91 > 0, 91, > 0, g > 0 and T > 
0, then 






47c"(H-5) 


Va'-i 


^ 4^.91, , . [t 1 ,1 « 

27ii'‘(H-5) W oi J 


Proof . Multiplying (1) by w* (w* is the complex conjugate of w) throughout, 
integrating the resultmg equation over the vertical range of z and utilizing (2) - (6), 
we get 


-a^){D^ -a^ w dz iQ -a^ - p*)Q*dz 
0 V 0 


•t91,o^ J(j)j -a^ |<|)*flfe-)-gj(£)^ -ct^)hz | -a^ -E-EL \h]dz 


-T\C, 


D^-a^-E- 


i;*dz-QT]i;\D‘^-a‘^-^^ldz = Q. ( 8 ) 


Subsequently, for convenience, we will omit the limits of integration from the integral 
sign and dz from the integrand. 
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Integrating the vanous teims of (8) by parts for an appropriate number of times and 
using the boundary conditions (7), we get 

J|(d^ +— +a^|M^p]+t9?,a^j([D(i)l^ +^3(^|<j)f ) 

+ Q\\{d^ -a%^^QTl\D^f +a^|^f]+7’j[i)Cf +a%f] 

-9?a'|[Def +a'|0|']+9?,aVJ|(j)|' +a^\Kt] 

-iRaV J|0f +^|;f + ^ 0. (9) 

Equating the real and imaginary parts of (9) to zero and canceling p, (^^O) throughout 
from the imagmary part, we have 

J|(o= +&|[d»|’ +a=|i»f ]+'t<R,a“j(jot,|' +a"|(t.f) 

+eJi(D* -a%f +QT\^f +a^|5f]+ +a’|?f] 
-9ia=jii)ef +a’|0f ]+,, ka'JW" *^i\Dh,f +o"|*,f ] 




( 10 ) 
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and 

+ Slto=|9|’-ij|C|’+2^J|5f =0. (11) 

Multiplying (1 1) by Pr and adding the resulting equation to (10), we have 
J|(z)a +^[f|z)>vp ] 

+ zRyi{D^" +a^|(t>f )+(2 J|(d' +QTj\D^^ 

+ 7’j’[z)Cf +a"|cf]-9la'j[£>e|' +a"|e|']=0. (12) 

Since, w, 0, 4>, hi and ^ vanish at z = 0 and r = 1, therefore Rayleigh-Ritz 
inequality’ yields 


(i) JlDu'f 

(ii) 

(iii) J|D<|)|^ 

(iv) J|D/j,|">7t’J|A,f 

(v) 


( 13 ) 


Further, utilizing Schwartz inequality, we have 
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>7t^J|vvj^ 


Consequently, 

This together with inequality (13(i)) gives 

J^jZ)^wj +la'^\Dwf +a‘*|ti'|" j ^ (tt^ +a^)^ J|w|^ 

Further, (2) implies that 

JkP = J>™’’ = J[(d= 

= ||(D'-c=>f+2;,j(D6|=+<,^|e|')+ipfJ|9f 

Smce, pr > 0, therefore it follows from (16) that 


and 




Also, emulating the derivation of inequality (14) and (15), we have 


[Using (13(i))] 


(14) 


(15) 


(16> 


(17) 


( 18 ) 
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J|(d^ + 2a’-iZ)0f +a^l0f j> (tc^ +«^)'l|ef (19) 

Combining inequalities (17) and (19), we have 

(20) 

Again 



><■'"'''* 56 ? 




fllfNr [Using (18) and (20)] 


>(71^ +a^)-< 1 + ^^ — ^ ► I- f0*(D^ Schwartz inequality] 

I v^'+^vj ' ' 




+a^ i 


Using inequality (15) in the first integral, inequalities (13) and (21) in (12) and 
utilizing the fact that pr ^ 0, we get 


(it’ +«’) JH’ +t«,<i’(it’ +o’)/w’ 

+e(it’ +o’)j(z«i/ +o%f]+T (it’ +<i’)j|Cf 
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^ 1 


()l= +a^J 


Eqn (1 1) upon using [13(i)] yields the following uiequalities 


fL.|2 _ Qo, r[^, |2 






JlCf (23) 


DKf (24) 


T%f>ie +a^)j|H'f -ea./Kr (25) 

Inequality (22) coupled with each of the inequalities (23)-(25) yields the 
following inequalities respectively : 


+r(,u4i-iW< «££LL i£r r 


(tt^ +a^f Jkr 


Ho') ,-iU|4,]= 


+T(»Ho')fl-AW<^ 


7c^ +a^ 


(n^+a^f 
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and 


2(71^ +<2(71^ 




(28) 


1 ^ 

Now, if 5 = min — , — , 1 , then depending on the value of 5, exactly one of the 

Oi J 

mequaUties (26)-(28) will imply that 





(29) 


Since the minimum value of 
follows from inequality (29) that 


-i ^ xxrr 


a‘ 


with respect to a is 


27%^ 


, It therefore 


27it^ 

4 


0 + 8 ) 



< IR = 


Inequality (30) implies that 

\p\ < (tc^ -I 


where 0 = 


( 491 ^ 

277t"(l-i-5)[~ 277t'‘(l-h5)J 


(30) 


( 31 ) 
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Further, it follows from inequality (29) that 

a' 


(32) 


Since the minimum value of 



with respect of is 47t^, therefore it 


follows from inequality (32) that 



(33) 


Combimng inequahties (31) & (33), we get 



471^1 + 5 ) 


-1 


This completes the proof of the theorem. 

Theorem 1, from the point of view of hydrodynamic stability theory, may be 
stated as : the complex growth rate p=Pr+ ip> of an arbitrary oscillatory perturbation 
of neutral or growing amplitude in magnetorotatory double-diffusive instability of 
Veronis type lies inside a semicircle in the right half of the prPr-p^sne whose centre is 

at the origin and whose radius is — -1 = — rP — -1 . This 

47C^(1+5) 47C^(1 + 5) 

result is uniformly valid for an initially top heavy as well as an initially bottom heavy 
configuration. 

Corollary 1 : If, (p, w, 0, ij), hi, Q, p = pi+ jp„ p, 0, is a nontrivial solution of 

eqn. (1) - (6) together with the boundary conditions (7) and 91 > 0, 9is > 0 and 

. 27 %* , ^ 

X < (l + 5), then pr < 0. 

491 , ^ ^ 


Proof : Follows from Theorem 1. 


I 
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Corollary 1 implies that oscillatory motions of growing amplitude are not allowed 
in magnetorotatory double-diffusive instability of Veroms type if the ini tial stability 

parameter X doesn't exceed the value (l-f-S). Further, this result is uniformly 

valid for an initially top heavy as well as an initially bottom heavy configuration. 
Remarks : The following remarks, now, deserve attention : 

277t^ 

(a) If 0 < 91 < 91s < — ^ — (l + S) and p, 0, then Corollary 1 implies that pr<0. 

2771^ 

(b) If 0 < 91 < — — (l +5)< 9is and p, ^ 0, even then Corollary 1 implies that 

/>r< 0. 

It is easy to see that 


(i) 6 = - for VTC, (ii) 5 = min 
a 




for hydromagnetic VTC, and 


(u) 


5 = min 


VO 


fox rotatory VTC. Consequently, the characterization 


theorem of Banerjee et al.^ and Gupta et al^ for VTC can easily be averred 
from (a). Further, in case of rotatory VTC and hydromagnetic VTC (a) 
yields more general results than those derived by Baneijee and co-authors*’^ 


The above remarks thus clearly and unequivocally establish the generality of the 
results derived herein. 

Theorem 2 (A semi-circle theorem for magnetoratatory STC) : 

If (p, w, 0, <|), hi, p = pr+ ip„ Pr ^0, p,^ 0, is a nontrivial solution of eqn. 
(1) - (6) together with the boundary conditions (7), and 91 < 0, 9lj < 0, then 
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Proof . Replacing 5R and 5Rj by -|5R1 and respectively in eqn. (1) and 
proceeding exactly as in Theorem 1, mutatis mutandis, we get the desired result. 

Corollary 2 : If (p, w, 0, (j), h„ l^,l)p=pr+ ip„ p, ^ 0, is nontrivial solution of 
eqn. (1) - (6) together with the boundary conditions (7) and 91 < 0, 91^ < 0 and 

i < (l + S), then pr < 0. 

4|9I| ^ ^ 


Proof : Follows from Theorem 2. 

The essential contents of Theorem 2 and Corollary 2 from the point of view of 
hydrodynamic stability are similar to those of Theorem 1 and Corollary 1 . However, 
now they pertain to magenotorotatory double-diffusive instability of Stem type. 
Further remarks similar to those after Corollary 1 hold here also. To be specific, 
Corollary 2 implies that . 


(a) If0<|91,|<|91|< ^^(U6)and/i,^0, then Pr<0. 

(b) If 0 < |9ij| < |9I| < — - — (l + 5] < |9i| and /?, A 0 even then pr < 0. 

4 ^ 


Further, it is easy to see that 


1 


(i) 5 = -forSTC,(ii) 5 = min -i,— 


1 1 




for hydromagnetic STC, and 


(iii) 8 = min 


'ill 

VCJ j 


for rotatory STC. Consequently, one can easily write down 


from Corollary 2 the characterization theorems for STC with or without the 
individual effects of a rotation and a magnetic field. 
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Abstract 

Lmear thermal stability of a horizontal fluid layer confined between two rigid boundaries and 
heated fi:om below is considered The temperature distnbution between the walls consists of the 
two parts a steady part and a time-dependent oscillatory part. Disturbances are assumed to be 
infinitesimal Only odd solutimis are considered. Floquet theory is used to study the stability of 
the solution Numerical results for the critical Rayleigh number are obtamed It is found that the 
disturbances are either synchronous with the primary temperature field or have half its 
firequency. 

(Keywords • modulation/ stability/ rayleigh number/ odd solution/thermal convection) 


Introduction 

This paper deals with the stability of a horizontal fluid layer confined between 
two ngid planes. The layer is heated from below and the temperature of the lower and 
die upper plates varies periodically with time. Study of convective flow is of great 
importance in geophysics and industrial operations. Considerable attention has been 
given to this problem during the last fifty years. Chandrasekhar' has given a 
comprehensive review of this stability problem with steady temperature gradient. 
Venezian^ who worked out the thermal analogue of Donnelly's ejqieriment^, was the 
first to consider the oscillatory temperature profile for modulating the wall 
temperatures, nevertheless a similar problem had been considered earlier by Gershuni 
and Zhukhovitskii'* for a temperature pro^le, obeying rectangular law. Later on many 
investigators studied this problem of thermal convection under various physical 
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conditions, considering different temperature profiles and boundary conditions. Some 
of them are : Rosenblat and Herbert^ Rosenblat and Tanaka®, Yih and Li’, Gresho and 
Sam*, Fmucane and Kell/, Roppo et al}'^ and Aniss et al}\ Recently author's 
group’’"'® has mvestigated the convective instability using sinusoidal function, saw- 
tooth function, step-function and day-night function for modulation of the wall 
temperatures. 

The objective of the present analysis is to determine the critical conditions under 
which convective instability can occur. The effect of temperature modulation on the 
stability limit of die system has also been investigated Only odd solutions have been 
considered. The results have their relevance with convective flows in the terrestnal 
atmosphere 


Formulation 

Consider a fluid layer of a viscous, incompressible fluid, confined between two 
parallel horizontal walls, one at r = -dll and the other at z=d/2. The walls are 
infinitely extended and rigid. Then using normal mode technique, the non- 
dimensionalized, linear governing equations, in the Boussinesq approximation are 


o (d’ -fl’)— = -P0+P(Z)’ 

dx 


( 1 ) 


0) 


fir \dz ) 


( 2 ) 


where ^ ~ perturbation quantities K=(m, v, w), 0 and p are respectively, the fluid 

velocity, temperature and pressure fields. T is the temperature in the conducting state 
while T is the time. P = vlk \s the Prandtl number, R = agATc^/vk is the Rayleigh 
number, a is the coefficient of volume expansion, g is the acceleration due to gravity, 
AT is the temperature difference between the walls, v is the kinematic viscosity, k is 
the thermal diffiisivity,'© is the non-dimensional frequency, and a is the horizontal 
wave number. 

The considered boundary conditions in the non-dimensional form are given by 
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r= 1-8 cos T at z = - — 

2 

(3a) 

= 8 COS X at z = — 

2 

(3b) 

where s is the amplitude of modulation. Then the temperature gradient 

be 

— is found to 
dz 

5^=-l+e RAFiz, X)} 

dz 

(4) 

where F(z, t) = ^2 = 

smh(A,/2) 

(5) 

The boundary conditions on w and 0 are 


w = Dw = 0 at z = ± — 

2 

(6) 

0 = 0 at z = ± — 

2 

(7) 

Method 

Expression (5) shows that F(z, i) is an even function of z. By carefully analyzing 
the eqns. (1) and (2) and the boundary conditions (6) and (7) we see that the proper 
solution of the eqns. (1) and (2) can be divided into two non-combining groups of 
even and odd solutions. Previous investigations'^’'* on thermal convection have shown 
that disturbances corresponding to even solutions are most unstable; however here we 
discuss the stability of the disturbances corresponding to odd solutions. 
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Now since 6 vanishes at z = ± , therefore it is expanded in a series of sin(2«7tz) 


Also w is written in a series of <])« so that 

(D^ - ay^„ = sin (2nnz) (8) 

where (}i„ = Dif'n =0 

Then the general solution (8) can be given by 

(|)„ = P„ sin haz + Q„ cos haz + sm (2n%z) ( 1 0) 

where P„ = (-lT - ?"^” -cosh(q/2) (11) 

sinha-a 

gn = -(-!)" ■sinh(o/2) (12) 

sinha-o 



(13) 

(14) 

(15) 


Now substitute (14) and (15) into the eqns. (1) and (2), and multiply by sin( 2 /M 7 tz) 
resulting equations are then integrated with respect to z in the interval -1 The 

I 2 2 } 
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outcome is a system of ordinary differential equations for the unknown coefficients 
A„(-z) and 


«£k 


-f 


(16) 


CO ciB„ 

__ fflt 

2 dx 


~[(2mny +Ra'^f[p„„ -b Re^„„e"}}4„ 

^ n=l 


(17) 


(«=1,2,3, ) 

The other coefficients, which cwcur in (16) and (17) are 


Pnn. = t!y2^A2)si^(2mnz)dz 

(18) 

=£j/^2-^^(j)„(2)sin(2/«Tcz)cfe 

(19) 

Lnm = £I/2'^'^<l>„(^)sin(2/«7rz)cfe 

(20) 

sin hiX/l) A(Xz)sm( 2 wn 2 )£fe 

(21) 


Here the values of the integrals (18)-(20) have been obtained in their closed forms, 
however (21) has been calculated numerically, using Simpson's (l/3)rd rule^®. Thus 

Pnm =^r«5„m +(-l)'”2»77tY*[2P„sinh(af/2) + g„{cosh(c3f/2)-4ay„sinh(fl/2)}] 

( 22 ) 

K„„ = -iyj(2«it)^5„„ +(-l)'”2»nty„[2(o^P„ +2og„)sinh(a/2) 

+a'^<2n{cosh(a/2) + 4£3ry„sinh(a/2)}] (23) 
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4m = |Y'(2«7t/5„m +(-ir2/K7iy«[2aV„ sinh(a/2) 

+ 'Q„ {2 - a^y„ )siiih(a / 2) + o" cosh(a / 2)}] (24) 

where 5„„ is the Kronecker delta. 

It is convenient for computational purpose to take m = 1,2,3 N i.e total IN 

equations and then rearrange them. For this, multiply the eqns. (16) by the inverse of 
the matrix (K„„-a^Pnm), and then introduce the notations 

xi-Ai,X2~Bi,X3=A2,X4 = B2 (25) 

Now combine the eqns. (16) and (17) to the form 

^^H„x,+H„X2+ +H^X, (26) 

at 

(/ = 1,2,3, 2Afandi: = 2N) 

where H,j (t) is the matnx of the coefficients in the eqns. (16) and (17). 

Analysis 

Since the coefficients H,j{x) of the eqns. (26) are either constant or periodic in t 
with period To= In/a, therefore the stability of the solution of (26) can be discussed on 
the basis of the Floquet theory^”. The solution in Floquet theory must be of the form 
eP^Pix), where P(x) is a periodic function of x with period Xq. The vanishing of the real 
part Pr of p gives the stability boundary. Here our aim is to determine the conditions 
under which pr=0. The disturbance is synchronoous with the unsteady part of the 
mean temperature field, if the imaginary part Pr=0. The disturbance is synchronous 
with the unsteady part of the mean temperature field, if the imagmary part p, of p is 
zero However if p,io is equal to ir or -it, then the disturbance is subharmonic, having 
frequency half that of the unsteady mean temperature field. 
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Let 


x„(x) = x,„(i) = co 1 [xi„(t),X 2 „(t), -Xjr„(T)] (27) 

(«=1,2,3, IN) 


be the solutions of (26) which satisfy the initial conditions 


x,„(0) = 5,„ (28) 

The solutions (27) with the conditions (28) form 2N linearly independent 
solutions of the eqn. (26). Once these solutions are found, one can get the vdues of 
x,„(2n) and then arrange them in the constant matrix 

C = [x,„(27r)]. (29) 

The eigenvalues X,;, X, 2 , X 3 , Xz, of the matrix C are also called the 

characteristic multipliers of the system (26) and the number defined by the 
relations 


Xr = exp (27tpr), r = 1,2,3, 2N (30) 

are the characteristic exponents. 

The values of the characteristic exponents determine the stability of the system. 
We assume that the pr are ordered so that 

Re(pi) ^ Re(p 2 ) ^ ^ Re(pi) (31) 

Then the system is stable if Re(pi) < 0, while Re(pi) =0 corresponds to one 
periodic solution and represents a stability boundary. This periodic disturbance is the 
only disturbance, which will manifest itself at marginally stability. 

To obtain the matrix C we have integrated the system (26) using Runge Kutta-Gill 

proceedu^e’^ The eigenvalues Xi, X 2 , X3, X^ of the matrix C are 

found with the help of Rutishauser method^'. 
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Results and Discussion 

By setting n, m equal to 1 and s = 0 (in the absence of modulation) in eqns. (16) 
and (17), we find the first approximation to Rayleigh number for the onset of 
convection. This minimum value of Rayleigh number at a particular wave number a is 
the value Rnait, at neutral stability This corresponds to sin(27tz), a trial fimction for 6. 
The corresponding value for R is 


R = - 


Uif+a^ 


-64oJt^ sinh^(fl/2)/|(45r^ +fir^)^(sinh a 


(32) 


which gives 

i2««/=17803.24 at 0 = 5.365. (33) 

This is m contrast to the exact value 17610.39, at the same value of the wave number. 
The second approximation to the Rayleigh number is found to be 17621.74 at 
a = 5.364, which is obtained by setting /», « = 1 and 2. Similarly the third 
approximation, obtained by putting n, m=l, 2, 3 is 17611.84 at o = 5.364 These values 
are same, as they should, to the Chandrasekhar's values*. By including more terms in 
the expansion of w and 0 one can achieve a higher degree of accuracy. With the help 
of the available packages one could carry out these calculations more accurately. 

The solution is obtained by solving the system (26) for X\, xt, x^ and X4. Now we 
calculate the modified value of Rc : (the critical Rayleigh number) with variation in 
other parameters and check the critical value of the wavenumber a, for the case when 
B ^ 0 Here the results are calculated for moderate values of e as we are interested only 
in the modulating effect of the oscillation, there seems to be no reason why this theory 
can not be applied for large values of the parameters. 

We obtain a relationship between the critical Rayleigh number Rc and the 
corresponding critical wavenumber in terms of the other parameters P, s and the 
dimensionless frequency o. The critical Rayleigh number is the minimum value of R 
as a fimction of the wavenumber a with variation in parameters P, s and ©. The 
corresponding value of a is known as the critical wavenumber Oc. To obtain the 
critical curve Rc versus e(Fig. 1) we proceed as follows: First fix P and to, then for 
some value of e, find that value of the wavenumber for which R is minimum This 
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minimum value of R and the corresponding value of the wavenumber are the critical 
Rayleigh number and the critical wavenumber respectively. Continue in this way we 
find Ro and the corresponding Oc for different values of e. Thus we get a curve Rc 
versus e, for fixed values of P and < 0 . The curve consists of two different curves, one 
corresponding to synchronous solution(S-curves) having same frequency as the 
applied temperature field and the other corresponding to subharmonic solution (H- 
curve) having frequency half that of the applied temperature field. 



Also we calculate the modified value of R,uut (R at o=5.364) with variation in 
other parameters and check the critical value of the wavenumber at neutral stability. 
Since R, is the minimiun value of R at some wavenumber (a^) while R„eut is the value 
of R at 0=5.364, therefore Rc is always less than R^t, but they coincide when B=0(in 
the absence of modulation). 

In Fig. 1 we have shown the variation of the critical Rayleigh number Rc, with s at 
<0=5.0 and P=0.13. The graph consists 5-curves and jy-curves alternatively, which are 
mtersecting at the cusps. Thus in the area above /f-curve there are also synchronous 
disturbances, but disturbances with half-frequency can be expected to be more 
unstable. Similarly above an 5-curve there are disturbances with half frequency but 
synchronous disturbances are more unstable. 
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It IS clear from the Fig 1 that initially in the synchronous region where 8 small, 
comection occurs at higher Rayleigh number than in the unmodulated case(8=0, 
Ra=\l()2l 74), so the effect of modulation on the stability of the system is stabilizing 
Also the cntical Rayleigh number increases as 8 mcreases so the modulatmg effect 
becomes more and more stabilizing. But on further increasing 8 we get alternatively 
/f-curve and S-curve and the cntical Rayleigh number Rc decreases. The value of 
decreases with mcreasing 8 and a stage is reached beyond which(e= 1.1) the value of 
R^ becomes less than 17621 74(unmodulated case), therefore convection occurs at 
lower Rayleigh number than in unmodulated case and the effect of modulation 
becomes destabilizing With further increasing 8, the modulation effect becomes more 
and more destabilizing as Rc decreases continuously. 

The results are in agreement with the findings of Yih and Le’. It also agrees with 
the results of Rosenblat and Tanaka^ and Bhadauria and Bhatia'** who found that the 
effect of the out of phase modulation of the temperature field is stabilizmg. 
Comparing the results of Rosenblat and Bhatia® and Bhadauria and Bhatia with the 
above results we see that their curves are smooth while here we have alternative S- 
curve and //-curve. Initially upto 8=0.28 we have a smooth curve and then alternative 
//-curve and .S-curve. As we go beyond a certain value of s we find a destabilizing 
effect, however in their cases the stabilization effect increases with 8. This difference 
may be because here the boundary conditions are slightly different from theirs. The 
above results also agree with the findings of Venezian^ and Bhatia and Bhadauna'^ 
who found that for out of phase modulation, initially the effect of unsteady part is 
stabilizmg but becomes destabilizmg at a particular P as co increases The results are 
in agreement with that of Yih and Li’ who found while studying the instability of 
unsteady flow that initially the effect of modulation is stabilizing and later on it 
becomes destabilized. They also found that the critical curve is composed of the two 
curves, one corresponding to synchronous solution and the other corresponds to 
subharmonic one. We also check the variation of the critical wavenumber Oc with 
respect to s, corresponding to the critical Rayleigh number Rc in Fig. 1. The critical 
value Qc appears to be discontinuous, as in Fig. 1, S'-curves and //-curves are not 
continued beyond their intersections. The existence of synchronous and subharmonic 
disturbances has already been indicated by Gresho and Sani®, Yih and Li’, Clever et 
al and Aniss et a/.” in their investigations. 

Fig. 2 depicts the value of the Rayleigh number (at neutral stability) with 
respect to 8 at a=5.36. Therefore the value of R^ut is slightly greater than the value of 
Rc in Fig 1. However the values of R in both the figures do not differ much, as the 
values of Uc corresponding to Fig. 1 are never very much different from 5.36. In Fig. 3 
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we consider the variation of with a at P = 0.73 and s = 0 1 The effect of 
modulation is found to be stabilizing, decreasing with increasing frequency a As a 
approaches infinity the effect of the unsteady part of the primary temperature field 
disappears altogether. These results are in agreement with that of Bhatia and 
Bhadauria'^ Bhadauria and Bhatia*'*, Rosenblat and Tanaka® and Venezian^. 



Fig 2- Variation of with e <0=5 0, F=0.73, a = 5 363808 



Fig. 3- Variation of with e co=0.1, P=0 73, a = 5 363808 


300 
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Abstract 

The problem of steady, laminar, free convection flow of a particulate suspension over an 
infinite, permeable, mclmed and isothermal flat plate m the presence of a transverse magnetic 
field and fluid heat absorption effects is studied An analytical solution is developed m the 
absence of particle loadmg and with constant particle phase density A parametric study 
depictmg the influence of the magnetic field and heat absorption effects is presented. The 
results for velocity, temperature and skin-friction co-efficients for both phases as well as the 
Nusselt number are illustrated graphically 

(Keywords hydromagnetic free convection/ particulate suspension/permeability) 

Introduction 

The study of free convection in porous media of an electrically conducting fluid 
past porous plate under the uifluence of a magnetic field has attracted many authors in 
view of Its applications to geophysics, astrophysics and controlled nuclear fusion and 
geothermal energy for power production. Moreau’ and other researchers have reported 
a large bulk of solutions and results for such flows. The study of dynamics of fluids 
with suspensions helps in imderstanding the scientific and engineering problems 
mvolved in pollution of city air, transport of suspended solid matenals m pipes, 
motion of aerosols in the upper atmosphere of the rotating earth etc. Cheng and his 
associates^ have extensively studied free convection and mixed convection about 
surfaces embedded in porous medium. In these analysis boundary layer 
approximations are used and Power law variations of wall temperature distribution are 
assumed Thangaraj^ studied on mixed convection about inclined surfeces in porous 
media with variable heat flux and Dandapat and Gupta'* worked on longwaves on a 
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layer of a visco-elastic fluid down an inclined plane Acharya et al.^ have studied 
effects of chemical and themial diffusion vntii Hall current on unsteady 
hydromagnetic flow near an infinite vertical porous plate. Dash and Rath® have 
studied applying explicit finite difference scheme for flow and heat transfer of an 
electrically conducting fluid Acharya et al ^ further studied a problem relating to heat 
and mass transfer over an accelerating surfaces. Recently Acharya et al} studied 
magnetic field effect on the free convection and mass transfer flow through porous 
medium. Mishra® have solved two-dimensional stagnation point flow of a dusty fluid 
near an oscillating plate. However, a few works are reported, so far as we are aware, 
about the hydro magnetic free convection of a particulate suspension from a 
permeable mclined plate for non-uniform particle-phase density. 

The object of this paper is to study the steady, laminar, free convection flow of a 
particulate suspension over an infinite, inclined and isothermal flat plate in the 
presence of a transverse magnetic field and fluid heat absorption 

Mathematical Analysis and Formulation of the Problem 

We consider the steady, laminar, hydro magnetic free convection flow of a fluid 
particulate suspension over an isothermal and permeable mfimtely long inclined flat 
plate Uniform fluid phase suction is imposed at the plate surface. Further an uniform 
magnetic field is applied normal to the flow direction. The fluid phase is assumed to 
be stagnant and is mamtained at a constant temperature Too, while particle phase is 
havmg an uniform density distribution and is affected by the gravitational force far 
away from the plate. 

The fluid is assumed to be Newtonian, electncally conducting and heat absorbing 
and retains constant properties except for the density in the buoyancy term. The 
induced magnetic and electric fields as well as the Hall effect of magnetohydro- 
dynamics are all assumed to be negligible and the magnetic Reynold's number, 
electrical displacement and convection are all assumed to be small Only the fluid 
phase will be affected by the presence of magnetic field since it is assumed to be 
electrically conducting However, the particle phase will not be affected by that field 
directly since it is assumed to be electrically non-conducting (insulator). But it will be 
influenced indirectly by the magnetic field due to the interphase drag mechanism 
between the phases. 

In the present study the particle phase is assumed to be made of non-uniformly 
distributed spherical particles all having one size and with a constant viscosity and 
diffusivity properties. The particle phase viscosity can be considered as a natural 
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consequence of the averaging process involved in representing a discrete system of 
particles as a continum’° " Also, the particle-phase viscous effects can be used to 
make particle-particle interaction and particle-wall interaction in relatively dense 
suspensions These effects have been investigated previously by many authors such as 
Tsuo and Gidaspow’^ and Gadiraju et al}^. Both the fluid and particle phases are 
taken as interacting continuously exchanging both momentum and heat transfer as 
discussed by Marble’"* The volume fraction of suspended particles is considered small 
compared to that of the fluid phase. 

The governing equations for the present problem are based on the conservation of 
mass, linear momentum, and energy for both the fluid and the particle phases. With 
above frame of references and assumptions (the plate is infinite in extent and the flow 
is steady), the physical vanables are functions of y only. Thus, the governing 
equations are given by 


dv 

dy 


( 1 ) 


d^u du 

dy dy 


dx 


■ ppN{u - « ) - pg cos (j) - ao^o « = 0 


( 2 ) 


c,N,{T,-T) + q^{T-TJ = 0 

dy^ dy 


(3) 


D 


d^Pp dippVj,) 




dy^ dy 


= 0 


(4) 


r ^ \ 


dy 


Pp 




“L 

dy 


-ppVp^ + Pp//(«-M )-p^gcos(|) = 0 

dy ^ 


(5) 


-Ppvp^ + ppNiv - Vp ) - p/.g sin (j> = 0 
dy\ dy ) dy 


( 6 ) 
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pp dy 



dy j 


PpCpVp 


tl^-p^c,Nr(Tp-T) = 0 

dy 


( 7 ) 


where 

Bo Magnetic induction, C: Fluid-phase specific heat, C/. Fluid-phase skin fnction 
co-efficient, D ■ Diffusion co-efficient, / • Any dependent variable, g . Acceleration 
due to gravity, Gr ■ Grashof number, H : Dimensionless gravitational acceleration, k 
Fluid-phase thermal conductivity, L : Characteristic length, M : Hartmann number, N ■ 
Interphase momentum transfer co-efficient, Np : Interphase heat transfer co-efficient, 
Ny Nusselt Number, : Fluid-phase pressure, Pr : Fluid-phase Prandtl number, qo . 
Dimensional heat absorption co-efficient, Qp : Dimensionless particle-phase density, 
Rv • Wall suctiqn velocity, S ; Dimensionless heat absorption coefficient. Sc. Inverse 
Schmidt number, T . Fluid-phase temperature, u . Fluid-phase x-component of 
velocity, U Fluid-phase dimensionless tangential velocity, v : Fluid-phase y- 
component of velocity, V . Fluid-phase dimensionless normal velocity, x,y . Cartesian 
co-ordmates, Y : Dunensionless normal distance, 

Greek Symbols 

a Velocity inverse Stokes number 
P Particle-phase to fluid Viscosity ratio 
Volume expansion co-efficient 
y Specific heats ratio 
s Temperature inverse Stokes number 
0 Fluid-phase dimensionless temperature 
K Particle loading 
}i Fluid-phase dynamic viscosity 
v Fluid-phase Kinematic viscosity 
p Fluid-phase density 
CTo Fluid-phase electrical conductivity 
<i» Tilt angle 

o Particle-phase wall slip co-efficient 
Op Particle phase 
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0,^ Plate wall 

()co Very large distance away from the plate surface (ambient condition) 
The boundary conditions are 


> 

1 

II 

o 


o 

II 

o 

w(oo) = 0 

ii 

T(oo) = T„ 

Pp(0) - pp^ 

pp(cc) = Pj=„ 

Up(0) - (0, -^y=0 

Up(co) = — ^ 

N 

Ii 

o 

Tj,(oc) = Too 


where coj is the dimensional particle-phase wall slip co-efficient. Smce the particle- 
phase may resemble a rarefied gas and undergoes slip at a boundary. Following Soo*’ 
and Chamkha'* we have specified the boundary conditions m the present problem 

The hydrostatic pressure gradient in eqn. (2) is approximated as 


dx 


= Pp„V«p„-p„gcos(|> 


( 9 ) 


where 


*Poo 


•gCOS({> 


( 10 ) 


This approximation can be easily obtained by evaluating eqns. (l)-{7) as y-^a>. 
Using Boussineq approximation to couple the fluid momentum eqn. to the 
temperature field and substituting eqns. (9) and (10), eqn. (2) can be written as 
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\id\ du 


dy^ " dy 


p^v——ppN{u-up) + p^ gcos^ 


+p„gp(7’-7’„)cos(j>-ao-Bo» = 0 
Introducing the dimensionless quantities 


Y = 


yGr* 
L ’ 


U = 


uL 

vgA 


V = 


vL 

vGr^ 


U - 

vGr^ 


jr 

^P=Vj> — r. 

vGr* 


T-T 
T -T ’ 

* W 00 


o JizL. 

T-T’ 


& = 


PP 

PpoO 






1 ’ 


k = Be2 
p« 


a = 


iVL^ 

vgA 


H 


_ SL^ 
v^Gr' 


vGr^ 


P=Vp/v 


2 _ 

i 


Af 


k 




\^p^p 

k, ’ 




Dp 


Gr = lMkzZkML 


( 11 ) 


( 12 ) 
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the eqns. (1,4) and (6) reduce to 


dY 

dY^ dY 

2^^^Qp^yQpV,^+aQ,iV-V,)+HGr^'^Q,sM = 0 
Now the dimensionless boundary conditions are 


0 

II 

1 

< 


C/(0) = 0, 

f/(oo) = 0 

e(o) = 1, 

e(Qo) = 0 

!2p(0) = Qpo, 

Qp{co) = 1 


Up{c^)=- 

0p(O) = 1, 

0p(qo) = 0, 


(13) 

(14) 

(15) 


(16) 


where Rv = L/vGr^'^ and a> = ajOr^'^IL are die dimensionless fluid phase suction 
velocity and the particle-phase slip co-efficient. Here all R^,QpQ and co are assumed 

to be constant. 

The fluid-phase skin friction co-efficient Cf, the particle-phase skin friction co- 
efficient and the Nusselt number Nu are defined, respectively as 
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du . 

= ryTi "rfF 


(17) 


dup j 


'/p 


1 




dY 


ir=o 


( 18 ) 


Nu=- 


hL dQ 


dY 


kGr 


1/4 


lr=o 


(19) 


Analytical Solution for Particle-Phase Density 
From eqn (13) with eqn. (16) we have 
V=-Rv 


( 20 ) 


Assuming Vp is constant throughout the domain of interest and substituting ecpi. 
(20) into eqn. (15) and rearranging we get 


Vp=—GAsm^-Ry = n 

a 


( 21 ) 


Obviously, for a vertical plate ((j)=0) both V and Vp wiU have the same wall suction 
velocity Rv. 

Substituting eqn. (21) into (14) we have 


Sc 


d^Qp 

dY^ 


dY 


( 22 ) 


The solution of eqn. (22) subject to the required boundary conditions can be 
shown to be 
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Qp=\+iQp^-\)e^^ (23) 

The solution (23) needs fiirther elaboration, which is very important A physically 
acceptable solution for Qp requires that it be always negative Eqn (23) bemg 
exponential, if ix is allowed to attain a positive value, it will grow towards infimty 
without limits; which is, obviously, not a valid solution. Physically, the requirement 
that Tt be negative means that the particle-phase normal velocity roust be negative i e. 
towards the plate, to ensure a contmuous flow of particles to compensate for the lost 
particles through the suction and to sustain the requirement of constant particle-phase 
density at the wall surface. If Vp is allowed to be positive i.e. away from the plate 
surface, all the particles will be cleaned away from the vicinity of the plate and the 
wall boundary conditions for Qp can not be met. The requirement of it to be always 
negative or at most zero was also required in the analj'tical and numerical solutions 
(see Ramadan and Chamkha*^’'®). However inorder to obtain a closed form solution 
we have studied the present problem m the absence of particle loading (i.e. A=0) and 
with constant particle phase density (i.e Qp = 1). In our future work a numerical 
solution is sought for to solve the coupled non-linear equations ansing out of due to 
particle phase loading {K^Q) and with variable particle-phase density. 

With above assumption and substituting eqns. (20), (21), (12) in eqns. (2), (3), (5) 
and (7) and on simplification we get 

+ -M^I7 + ecos<|) = 0 (24) 

dY^ dY ^ 


+ 50 = 0 


dY^ 


dY 


(25) 


p + a(i7 - f/p ) -77 cos (j) = 0 

dY^ dY 


BPr 


nfZ£.+8(0-e^) = O 
dY 


(26) 


(27) 



450 


G.C. DASH and N. NAYAK 


Solution 


e = e-'” 2 ’' 


U = [e-^r _g«4r 1 

- 2 n 1 J'2 L J 




If e”^ 

U^=C^e’^ --cosif^—^ A, 

a ^2 ® 




2 71 a 

/K 4 7774 

2 P ^ P 


0 p=e'^+- 


ml +ff»2 ^Pfp ~^P^P 


(28) 

(29) 

(30) 

( 31 ) 


where the constants, C2, mi>Q, 7774 < 0 , TTie < 0 , 7773 < 0 are given in the Appendix. 

Results and Discussion 

Following numerical values were assigned for various physical parameters during 
computations. 


G7=10,000, Pr=0n, Prr-2, J?v=0.5, /f=l, a=l, p=0.5, e= 1, (|)=15‘, y=l, 0=2 


From the Fig. 1 , it is observed that as the square of the Hartmaim number 
increases both fluid phase (solid lines) and particle phase velocity (broken lines) 
decrease But it is interesting to note that particle phase velocity profiles assume 
negative values for higher values of Hartmann number. It is further observed that in 
the absence of magnetic field the particle phase velocity becomes negative near the 
plate wall as well as far away from the wall. Consequently particle phase velocity 
leads to flow reversal effect due to increase in magnetic field strength. Therefore it 
may be concluded that increasing the Hartmann number leads to damping of the fluid 
and the particle phase velocity profiles. This is due to the application of a transverse 
magnetic field normal to the flow direction which imposes a resistive force (Lorentz 
force) similar to the drag force which tends to resist the fluid flow and thus reducing 
the velocity. 
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u, u. 



Fig. 1- Effect of Hartmann number (5=0); U : Fluid velocity profiles {Solid lines); 

Up Particle phase velocity profiles (Broken lines) 

It is also recorded from the numerical values of the temperature profiles that for 
different values of Hartmann number temperature profiles for both fluid and particle 
phases coincide, hence graphs are not presented here. Hence it is concluded that 
Hartmann number has no effect on the temperature profiles of the fluid and particle 
phases. Thus the Hartmann number has no effect on the Nusselt number. 

The effect of heat absorption co-efficient S, are depicted in Fig. (2 & 3). It is 
observed that flie velocity and temperature profiles for both fluid (solid lines) and 
particle phase (broken lines) decrease when S assumes the values 0,-1 and -3. Hence 
it is concluded that the effect of heat absorption is to damp the flow and heat transfer 
phenomena of both phases and to reduce the thickness of both the velocity and the 
thermal layers. This result is realistic one since, due to heat absorption, less energy is 
used to support the buoyancy effects while the considerable amount of energy is 
absorbed by the fluid. The above phenomena jure clearly established as U, Up and 0 
decrease as j5| increases. 
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Fig 2~ Effect of heat absorption (M=0), U * Fluid velocity profiles (Solid lines). 
Up . Particle phase velocity profiles (Broken lines) 

o. o. 
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The present study couldnot be undertaken in the presence of heat generation co- 
jfficient (S > 0) as the solution of eqn (26) does not exist 

The effects of magnetic field (solid line) and the heat absorption (broken line) on 
he skm-ffiction of fluid and particle phases. C/and C/p, are shown in Fig 4. From the 
5 qn (18) it is seen that Cjp is identicdly zero in the absence of particle loading, (K=0) 
Ihe decrease of C/ with the increase of magnetic field and heat absorption is noticed 
rhis is due to the dampmg effect produced by magnetic field and heat absoiption 
larameter on the flow. 


Cf 



Fig. 4- Flmd-skin friction coefficient; effect of Hartmann number (5=0) (Solid line), 
effect of heat absorption (M-Q) (Broken line) 


Lastly, the effect of heat absorption on the Nusseit number is shown in Fig. 5. The 
'Jusselt number increases with the increase of heat absorption parameter. This may be 
ittributed to the feict that as the heat absorpticm increases the plate tends to supply 
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more energy through its wall to maintain the co n sta nt temperature of the wall, and 
hence the Nusselt number. 

Nu 



Fig. 5- Nusselt number (Effect of heat absorption (A/=0)) 

Conclusion 

The mathematical modelling for free convection flow of a particulate suspension 
over an infinite, inclined, permeable and isothermal plate in the presence of both 
magnetic field and fluid heat absorption was under taken. The present study accounted 
for both particle viscous and diffiisive effects. The general effects of the magnetic 
field and the heat absorbing fluid were found to dmp the flow and the thermal 
profiles. Moreover, the Nusselt number was found to increase with the heat absorption 
effects but the Hartmaim number has a little effecTbn the Nusselt number. 

Appendix 


1712 =PrRv+4Pr^Iiv-4PrS 
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Abstract 

In this bnef write up, we indicate without proofs the interesting results connectmg the troika of 
the topics in the title after giving their defimtions and some examples This article is by no 
means exhaustive and it reflects the interests of the author 

Order 

The field of rational numbers has the binary relation of order besides the two 
algebraic operations (addition and multiplication). This ordering is compatible with 
these algebraic operations. The significance of order came to the fore only in the 19th 
century with Cauchy and Dedekind independently developing the completion of the 
rational number field to real nximbers, wluch lent rigour in real analysis, in different 
ways using order. Hilbert contributed to the resurgence of interest on order when he 
developed an axiomatic approach to geometry. However, an impetus given by Artin 
made the study of ordered fields a topic of great mterest and utility. 

We begin with the notion of an ordering P of a field K. This ordering P is subset 
of K havmg the following properties. P + P c. P, P.P c, P r\ {-P} =. {0} and 
Pu{-P} = K It is clear that P contains (finite) sums of squares of elements of K. 
Generalizing this notion we define a preordering P in a commutative ring R with 
identity, as a proper subset P of P such that c P, and P + P ^P, P.P c P. The 
requirement of P as a proper subset of R may be strengthened to -1 g P. 

From now on we shall restrict ourselves to the case when P is a field K From our 
definition of a (pre) ordering P in AT, we can immediately infer that the characteristic 
of K is zero (for otherwise, the characteristic is p 0 and -1 = 1 + 1 + ...+ !, (p-1) 
times which shows -1 e P, a contradiction). We also find that -1 is not a sum of 
squares. Given this P which is called a positive cone in case P is an ordering, we 
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introduce an order < by setting a ^ b in case b — a € /*. It is easy to verify that with the 
above stipulation of order, K is an ordered field. This is in the sense that the field K 
has a linear order (so that any two distinct elements are comparable) which is 
compatible with the addition and multiplication by positive elements ofK. 

Some examples of ordered fields. 

1 The field of rational numbers Q and its completion with respect to the 
topology induced by the metric arising out of ordinary absolute value (also 
known as the Euclidean topology), namely, the field of real numbers R, have 
only one order 

2. The field Q{t) of rational functions over Q in one variable t has uncountably 
many orders. This is obtained by imbedding Q{t) in the *real field R with t 
being mapped onto some transcendental element in R over Q. 

3. Given the field of rational numbers g ©r the field of real numbers R, any 
transcendental extension <2(0 {Rit)} can be ordered by choosing as the set of 
positive elements those which are of the form o’.t7(t).{g(0}~’ where a is a 
positive element m Q{R),f{t) and g(t) are polynomials with constant term 1. 
(See Jacobson [1], pp. 272.) 

Given a preorder f on a field K (or even in a rmg R) we can extend it to a larger 
preorder, T' containing T. For this, suppose o is a non zero element K with -a not m 
T. Then T' = T[d\ which is the sub ring of K generated by T and a satisfies the 
conditions of a preorder In fact a necessary and sufficient condition for this T'to be a 
preorder is that as -T. Now the set of all pre-orderings P under set inclusion forms an 
inductive set Using Zorn’s lemma we find that any maximal element in P is indeed an 
ordermg of K (or R). 

ArtSi defined a field to be formally real if -1 is not ejqiressible as a sum of 
squares in K. We have the following ; 

Theorem 1 : Let ilT be a field. Ihen K is formally real if and only if AT has an ordermg. 

For the proof, observe that if K is formally real, then taking T as the suifi of 
squares of elements of AT, we get T as a preorder in K and this can be extended to an 
order as ejqilained in the previous paragraph. For the converse, note that if K has an 
ordering then -1 is not a sum of squares. 

Artin was led to the study of formally real fields while solving the seventeenth 
problem of Hilbert's celebrated problems. Hilbert^ and Landau^ had shown that if a is 
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an algebraic integer which is positive in every real embedding of Q{a) in R then a is 
expressible as a sum of squares in Q(a.). In fact if F is a finite extension of th^ field of 
rational numbers and {o, 1 1 = 1 , 2,...r} are the different embeddings of F in the field of 
all real algebraic numbers, then any non zero element 7 € F is expressible as a sum of 
squares in F if and only if o, ( 7 ) < 0 for all /. 


Hilbert had shown earlier that if a rational function (p(x) e Q(x) {(p(x,y) e Qix,y)} 
assumes nonnegative values whenever the denominator of (p(x) {or (p(x,y)} does not 
vanish, .t|i^ <p(x) {^(x.y)} can be expressed as a sum of 2(4) squares 
. The 17th problem is die following : 


1-1 


Z/X^f 'ZfMy) 




If (p(xi, X 2 ,.. x„)mK = Q(xi, X 2 ,...Xn) is positive semi definite whenever it is defined 
over g, is (p a sum of squares in AT? The quantitative aspect of this was conjectured by 
J.Ax in 1968 and proved by A. pfister'* (positive semi definite means (p(oi, 02 ,. ..a„) > 0 
where a\ e Q. 


Artm solved the problem (the qualitative version) in the affirmative and for this he 
mtioduced and with Schreier characterized the real closed fields. A real closed field is 
a formally real field, which does not allow of any formally real proper algebraic 
extensions. These real closed fields have the property that any polynomial of odd 
degree (with coefficients from the real closed field) has a root in it. Moreover the 
adjunction of square root of -1 gives the algebraic closure of real closed field. In a 
real closed field, the positive cone is given by the sums of squares of its elements. 
Artm^ proved the following 

A 

Theorem 2 Let (K,P) be an ordered field with K as its real closure (which exists and 
IS unique up to AT-isomorphism). Denote by K„ = K{y\, X2,...x„) the rational function 
field m n variables over K Let f e K„\k such thaty(“) <^ 2 , ■"») ^ 0 in AT for all a 

i 

e (oi, 02, ..Op) e IC such that/xi, X2,...x„) = ^ bj^ for a suitable integer;. 

1=1 


In case P is the only order of AT, then b, is a sum of squares in K for each ?. (For 
instance, K is an Euclidean field or a Pythagorean field). 

We will not give a proof of the above theorem. One may find a proof in Jacobson* 
or Lorenz®. Artin used Sturm's theorem in counting the real zeros of a polynomial over 
the real field in a given interval and the theory of formally real fields. 
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Valuations 

The origins of the theory of valuations may be traced to Kummer. It came to 
prominence with the work of Hensel when he investigated solutions of congruences 
modulo different powers of a prime p. Other major contributors are Rychhk, 
Ostrowski, Krull, Hasse, etc This theory also gained, from an algebraic approach to 
geometry, especially when one considers the behavior (or properties of an algebraic 
variety) at a point Thus the notion of a place from geometry is equivalent to 
valuations and valuation ring as we see below. We begin with the definition of general 
valuation as given by Krull Other (equivalent) definitions are due to Artm^ for rank 1 
valuations and to Nagata* for general valuations. 

Definition 1 . A valuation v on a field is a homomorphism of the multiplicative 
group K* of nonzero elements of /T to an ordered abelian group F with the additional 
property 


v{a±b) > Min {v(a), v(6)}. 

By settmg v(0) = oo, where oo > y for all y in F and oo having the usual properties 
associated with the symbol, we extend v to K. 

The following properties are easy to verify : 

(a) 1{D= {a e k\ v(a) > 0} then Z) is a subring of K. 

(b) The set M = {a e I v(a) >0} is a maximal ideal of D. 

(c) All elements of D not in M are units in D. 

(d) If a,bmK with v(a) < v(b), then v{a+b) = v(o). 

The ring D is called the valuation ring of K with respect to the valuation v. It is a 
local ring with M as its unique maximal ideal. The field DIM is called the residue class 
field of the valuation v. Here are some examples of fields with valuations. 

(i) Every field K has the trivial valuation namely v(a) = 0 for all a 0 in AT and 
v(0) = 00 . 

(ii) In the field of rational numbers Q, for each prime p, we define Vp(r) = « if /?" 
divides r exactly 

This means writing r = in the lowest terms so that (a,b) = \, p” divides a exactly 
and (p,b) = 1, arp" divides b exactly and (p,a) = 1 (in case «<0). 
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(hi) The field of rational numbers Q, real numbers R and complex numbers C 
have the usual absolute value (;.e. \a\= a \ia>0 and -o if a < 0 for all a e 
Q or R and \a\ = positive square root of N(a) = norm of a = if a e C 

{a = X + j.y). This associates with each element of (Q, i? or C) an element of 
the set of positive real numbers which satisfies the following properties • 
\a\ ^ 0 and equality holds if and only if a = 0, \a.b\ = |a|.|6| and \adi)\ < \a\ + 
\b\. This fimction arising out of the absolute value is known as the 
Archimedean valuatidn Note that this valuation satisfies a weaker inequality 
than the one given in the definition and takes the values in the multiplicative 
group of positive real numbers and zero added to this group. 

(iv) In the field of rational functions K = k{x) in one variable x over a field k, 
each irreducible polynomial p(x) in ^[x] gives a valuation (this is given 
exactly as in example (iii) with the irreducible polynomial playing the role 
of a prime and a and b are polynomials). Besides these, there is the valuation 
given by x"\ We now define the notion of a place. 

Definition 2 • A place on a field A" is a homomorphism (p of a subring of K mto a 
field k satisfymg the following: 

(a) if a e K, a i then e Kq, and q)(o'*) = 0. 

(b) there exists some x in for which q>(x) 0. 

This notion comes from algebraic geometry. It can be shown that if is a subring 
of a field K and .4 is a proper ideal of R, then there exists a place (p of K such that 3 
R and^ = {a e K \ (p(a) = 0}. 

Examples of places : 

1 Consider an irreducible element / m a unique factorization domain R. Then 
the ideal A generated by /is a prime ideal so that R/A is an integral domain. 
Let K be the field of quotients of R and k the field of quotients of R/A. Then 
the canonical homomorphism (p of y? to R/A can be extended to the subring 
Rf consisting of all those fractions a/b with a, b in R and b e A. This 
extended homomorphism is indeed a place. 

2. If A is a Dedekind domain, and Af is a maximal ideal of A then the 
homomorphism of A onto A/M can bt extended to the local ring Am = {ct/b, 
ab e A, b^M} e AT the field of frac’;ions of A. This extension is again a 
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place (For more information about places see Zariski-Samuel's 
Commutative Algebra. (Vol. II, Chapter VI). 

We have the notion of a valuation rmg, which is defined as follows. 

Definition 3 : A valuation ring i? of a field AT is a subring of K such that of x and x'' 
(x 0) in K, at least one belongs to R. 

Examples of valuation rings. 

1. Consider the field of rational numbers Q with the valuation Vp given a prime 
p. Then 

Zp= {r = a/b e Q\p does not divide 6} is a valuation ring. 

2. In the field K = k{x) of rational fimctions in x over k, let p(x) e k[x] be an 
irreducible polynomial of degree > 0. Then Rp -^{f/g e AT |/, g e k[x] and 
/>(x) € ^[x] does not divide g(x)} is a valuation ring. 

3 . Any finite field F has only the trivial valuation. This is clear since any non 
zero element in F is a root of unity. 

We are ready to show that the three notiois valuation, place and valuation ring are 
equivalent. 

Let (^,v,r) be a field with a valuation with value group F (with oo added to it). Let 
i? = (a e A’ I v(a) > 0}. Then ii is a ring with identity. If a e K does not belong to R, 
then v(a) < 0 so that -v(o) = v(a”*) > 0 and therefore o"’ e R. Thus F is a valuation 
ring. 

Now given a valuation ring /? of a field AT, consider the ideal M consisting of those 
a m a: which are in F while a'‘ « R. This M is a maximal ideal ofR. Setting k = R/M 
and defining a homomorphism 9 from AT to i by (p(a) = a as the canonical image of 
a, in i in case o e J? and 9 ( 0 ) = 00 if o ^ we see that 9 is a place on K. 

Finally, given a place 9 of AT whose image is in the field k, we shall produce an 
ordered abelian group F and a valuation v on AT with value group F. For this purpose, 
consider the multiplicative group of non-zero elements of AT and the set 17 = 9 '’()t*) 
the pre image of the non-zero elements of it. 

Denote the quotient group t^UhyV fixe natural homomorphism from AT* to F by 
V. We order this F as follows. The positive part F+ of F is the image of the setAf under 
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the map v where M is the set of those elements of K* which are mapped onto 0 by cp. 
(We shall use the additive notation for the group operation mduced by the 
multiplication in K). As this set M is closed under multiplication the set r+ is closed 
under addition, that is to say r+ + r+cr+. 

Let X € K*\M. Then v(x) is not in r+ If a € F there exists an a in such that v(a) 
= a Should a belong to r+, then o is in M so that a~^ g M and v(a("') = -v(a) = -a g 
r+. Thus r+ n (-r+) = {O}. in case a e r+ and a then a g t/ and a g M so that 
a”’ e M and therefore -a e r+ fiilfills the conditions for the positive part of a totally 
ordered group. 

To establish that v is a valuation we need to show that v(a±b) > Min {v(a), v(b)} . 

Suppose v(a) > v(b). Then v (1 + a/b) = v(a+b) - v(b) we observe that a/b and 
therefore also (a+a/b), belong to the pre image of r+ u {oo} which shows that 
v{a+c b) = v(p+b) - v(b). Thus v{a+b) t Min {v(ar), v(h)}. From the way we have 
obtamed v and F it is obvious that v is onto. Thus {K*, v, F) is a field with a valuation. 

As mentioned earlier, the fields of real and complex numbers have limiting 
processes due to the presence of the archimedean valuation arising out of the usual 
absolute value and the usual metric determined by this valuation In a similar manner 
the p-adic valuation for different primes /j in the field of rational numbers gives rise to 
ultra metrics (as these metrics satisfy a stronger inequality) which are different for 
distinct primes. Each of these /?-adic topologies results in different limiting processes. 
Yet, the exceptional position of the topological fields of real and complex numbers is 
not an accident but a necessity arising out of general considerations. In this 
connection, the following theorem of Pontrajagin is important. 

Theorem 3 ; Let AT be a locally compact field, which is connected and satisfies the 
second axiom of countability. Then K is isomorphic to one of the following fields : the 
field of real numbers, the field of complex numbers or the skew field (in case K is not 
commutative) of real quaternion. 


Quadratic Forms 

The history of quadratic forms can be traced to antiquity. For instance, the Greeks 
knew the Pythagorean triples. The Renaissance came about with the theorem of 
Lagrange, which says that every natural number (and therefore also every positive 
rational number) can be expressed as a sum of at most four squares and that there are 
numbers which require four squares. During the 19th and early 20th century a good 
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deal of development took place in the arithmetic theory of quadratic forms. For 
mstance, in the first quarter of 20th century, C.L. Siegel developed the arithmetic 
theory of quadratic forms, which carries over from the field of rational numbers to 
algebraic number field and much more Among several applications of his results the 
thing that is relevant for our discussion is the following. IfX^) = 2(Or cos 
sing r&) where r runs from 1 to n, a„ b, are real numbers, is non-negative definite (that 
is to say /Q) > 0 for all values of 9), then y(9) has a representation as a sum of 
squares. In fact, it can be expressed as a sum of eight squares with rational coefficients 
and trigonometric polynomials in 9/2. In other words, 

M) = = l{<cos(r .9/2)-h< sin (r.S/2)} 

in which c*, are rational numbers, r -aries from 1 to « and s varies fi'om 1 to 8. 

The above result may be viewed as an extension to trigonometric polynomials of 
the followmg theorem of Hilbert^ and Landau^ and mentioned in the section on Order. 

Theorem 4 : Let A!’ be a field of characteristic different fi'om 2. Then an element a in AT 
is a sum of squares if and only if a is totally positive. Here by totally positive we 
understand that if K is finite extension of the field of rational numbers then a in AT is 
positive in every real embedding of AT and in case AT is a finite field then every a is 
totally positive 

An «-ary quadratic form over a commutative ring i? or a field AT is a homogeneous 
polynomial of degree 2 in «-variables. (We assume that 2 is not the characteristic of 
K) It may be expressed 3sJ{xx,..x„) = la,jx;cj, I <i,j^n. Since the characteristic of AT 
IS not 2, we may replace a,j by {a,j + ay)/2 in which case the matrix representing / 
namely the matrix of coefficients of /is symmetric. 

Given an w-ary quadratic form g, we define a quadratic space V over AT as an n- 
dimensional vector space over k with the map (also denoted by) g fi'om V to K 
satisfying the following properties 

(a) g(a.v) = a^.g(v) for ain ATand vin F. 

(b) The map : F X F -> a: given by A(v,,V 2 ) = - {g(v, +V 2 ) -g(vi) -^(vz)} is K- 

bilinear. 2 

Let (Fi, gi) and (F 2 , gz) be two quadratic spaces over AT. We say these are 
isometric if there exists a AT-linear isomorphism T of Fj to Fz such that gi(v) = g 2 (Tv) 
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for all V m Fi It is not difficult to see that the isometry classes of n-dimensional 
quadratic spaces over A" are in 1-1 correspondence with the equivalence classes of «- 
ary quadratic forms over K. (By equivalence of two n-ary quadratic forms gi and g 2 
we mean the following : there exist a nonsingular linear transformation T in GHn,K) 
such that gx(X) = g 2 iTX) whereXis the row vector (Xi, X 2 ,..X„) of the «-variables). 

In the late 1930's Witt developed the algebraic theory of quadratic forms by 
introducmg the orthogonal sum and Kronecker (or tensor) product in the collection of 
eqmvalence classes of quadratic spaces over K. Thus, if (F, gi) and (F 2 , ^ 2 ) are two 
quadratic spaces over K of dimensions «i and «2 respectively, then the quadratic space 
(F, g) of dimension («i + n-^ is given by F= F] © Fi and g(v) g(vi+ V 2 ) = g(vi) + g 2 (v- 
2 ) where v = vi + V 2 . If A, is the symmetric matrix representing g, for i = 1, 2, then Ae 

. It can be seen that the orthogonal sum depends only on the 

eqmvalence classes of the summands. 

If, instead of the direct sum of Fj and F 2 , we take the tensor product of these two 
spaces over K namely If ® If and take the quadratic form g detennined by the matrix 
Ai 0 A 2 , we get (F,g)as a quadratic space of dimension «i, «2 and this dso depends 
only on the equivalence class. We remark that if g, is a diagonal form (so that the 
matrix A, is a diagonal matrix) then g = gi ® g 2 is given by 





Thus if 

II 

S2 = {i>xA>A,) 

then 

g = gl -L g2 = (^i,. 


and 

gi ® g2 =(...a,b,... 

> 1 ^ /■ < «i., 1 <J < «2 


The set of equivalence classes of quadratic forms over K forms a cancellative 
semi ring under the above operations of orthogonal sum and Kronecker product. Here 
we need Witt’s cancellation Aeorem namely. If g, gi, g 2 are three quadratic forms over 
K such that g © gi s g 2 (over K). Here the empty form g = 0 of dimension 0 acts as the 
zero element and the form 1 = (1> of dimension 1 as the unit element. The semi-ring 


matrix of g is 


A, 0 
0 A, 
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of equivalence classes of quadratic forms can be extended as a ring and this is the Witt 
ring of i:, and denoted by WiK). (For details one may refer to T.Y. Lam. 191). 

Examples ofWitt-nngs. 

(1) Let /T = i? be the field of real number. Then the only anisotropic forms over 
K are’O, n <1) and « <-l> where n is any natural number. Thus W(R) = Z. 
Here a quadratic form is said to be anisotropic if it does not represent 0 non- 
tnviaUy. 

(2) Let K be any field of finite characteristic p. Say K = ZJpZ where p is an odd 
prime. Then we have W (ZIpZ) = {Z/IZ) © {Z/IZ) in case p = l(mod 4) 

and W {ZIpZ) = (Z/4Z) m case p s (mod 4). 

(3) If K is any algebraically closed field, or K is quadratically closed with 
charactenstic different from 2 so that every element is a square, we have 
W{K) = {Z/2Z). 

(4) From (2) above, using Hensel's lemma or otherwise one can show that 

= {ZnZt if/> 3 1 (mod 4) 
and W{Qp) = (Z/4Z)^ if p = 3 (mod 4). 

Relationship between Order and Valuation 

In this section we shall investigate the valuations that arise naturally on a field 
with ordermgs. We will therefore be interested in orderable fields. Given a formally 
real field K, there is a natural valuation on K whose residue class field is also a 
formally real field This natural valuation is obtained as follows ; 

As .ST is formally real, it contains a sub-field isomorphic to the field Q of rational 
numbers and also it is orderable. Let P denote the positive cone of some ordering of 
K. 


Set A{Q,P) - R - {a in 1 for some r ^ Q, r ± a e P} and 
= {a in K\r±a^P for all r > 0 in (2} 
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In can be verified that A(Q,P) is a ring and M{Q,P) is its unique maximal ideal. 
This ring is the convex hull of Q 'mK with respect to the chosen ordering P, that is to 
say if n < c < 6 and a,b e R then c e R (because r + c = r + a + (c-a) belongs to P as 
(c-a) € P and r±aeP and r-c = r-b + b- cis'mP, since r±beP and {b-c) e 
P Further, the maximal ideal M is also convex in R. For this suppose 0 < a < b, b e 
M Then we have 0 <b~^ < d'^. Since b e M, b~^ £ R and as i? is convex, a~^ ^ R so 
that a belongs to M. The residue class field R/Mhas the induced order with respect to 
which the positive cone is the canonical image of P i? in R/M Thus the valuation v 
of K, whose valuation ring is R and residue class field which is an orderable field, 
is a real valuation arising out of the order P vciK. We call this valuation nng A{Q,P) - 
R, the canomcal valuation ring of P and the associated valuation v (denoted by Vp to 
signify that it anses from P) as the canonical valuation. 

This canonical valuation Vp is compatible with the ordering P of P in the sense 
tliat if Q<a<b with respect to the order P, tjien v(a) t v{b). 

We have the following : 

Theorem 4 • Let P be a given ordering of the field K. Then the family F of all 
valuation rmgs in K which are compatible with the ordering P forms a chain under set 
inclusion. The smallest member is the convex hull of Q with respect to P in AT namely, 
the ring A(Q,P)- 

In fact, F consists of all sub rings of K which contain .4 (g.P). 

The residue class field RIM of the canonical valuation Vp of K for the given 
ordering P is in fact a real valuation. The field R/M with the push down order p (say) 
is an archimedean ordered field (because for any x e P since -r<x<r with respect to 
the ordering P, for some r in P, by pushing down to R/M we have -r < x < r with 
respect to the order P ). Consequently {R/M, P) has & unique order imbedding into 
the field of real numbers with its usual ordering. 

The above notion of canonical valuation Vp of an ordered field (K,P) is implicit in 
the paper of Artin-Schreier and Baer and explicit in the paper of Ktu11'°. 

Theorem 5 : Let (K,P) be an ordered field, P its canonical valuation ring for P and v 
be the canonical valuation. Then die following are equivalent. 

(1) There exists a real valuation ring S contained in the real closure K of {K,P) 
such that Sr\K = R 
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(2) There existe an ordering P of the quotient field K of the valuation ring 
S^- A(Q, P) of K . This valuation ring is the Henselization of {R,v) and 
Pr^K = P. 

Given a real field K, we can consider the set Xk of all orderings in K. This set is 
clearly non-empty as the pre-ordering (which may turn out to be an ordering) arising 
out of the sums of squares can be expanded to a positive cone of K. We can make Xk 
mto a topological space by using the "Harrison Sets" 

H(a) = {P e X\a €. Xk} for any a e K, a * 0, 


as a sub basis. With this topology Xk becomes a compact Hausdorff space and this is 
also totally disconnected. This topological space of orderings of a field and its ring 
theoretic analogue are used in the study of real algebraic geometry. 


Suppose now we have an orderable field K that has a valuation v. We may 
consider of the orderings of K, which are compatible with the valuation v. The 
natural question is what possible relationships can exist between and the value 
group r of V? 


Note that if a in X is a square, say a = b^, then v{a) = 2v{b) so that v(o)e2; 
(which by abuse of notation we write as F^). We then get a mapping v from 



to 



. If we take v to be onto, the v has a cross section. Thus there 


exists a map |i from 




such that v.p = identity. 


The theorem we state below without proof gives a relationship between Xk 
and r. 


Theorem 6 . Let X be an orderable field with valuation v and value group F. Then 
there is a 1-1 correspondence between the space Xk of orderings of K that are 


compatible wWi the valuation v and the set x Xt where ^j) = 


Horn 


and Xk is the space of all orderings of the residue class field k of AT with 


respect to the valuation v. 
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An immediate corollary of the above is the following result. 

If the value group F is the 2-divisible so that F = F^, then every ordering of k has a 
unique lifting to an order of K which is conqiatible with v. 

Without the above restriction on F, we see that each ordenng of k has as many 
lifting as the index of F^ in F to orderings in K which are compatible with v 

Relationship between valuations and quadratic forms : 

We begin this section with the well-known local global principle for quadratic 
forms known also as Hasse-Minkowski theorem. 

Theorem 1. Let A" be a global field of characteristic different from 2. Lety(xi,... ,x„) be 
an nensional quadratic form over K. Then / represents zero non-trivially over K if 
and only if, / represents zero nontrivially over all completions Kp where p is non- 
archimedean or real. 

The above theorem, also known as Hasse-Minkowski's local global principle can 
be stated as follows : two quadratic forms and gz over a global field K are 
equivalent if and only if, they are equivalent over for all prime ideals p including 
the infinite primes. Here denotes the completion of K with respect to the valuation 
arising out of the prime p 

We will now record some results connecting complete discrete valuation rings and 
their Witt rings. Let R be a discrete valuation ring with maximal ideed M and residue 
class field R/K = k. We shall assume that 2 is a unit in R. Then in case R is complete 
under this discrete valuation. R is 2-Henselian. This means that o in /? is a square if its 
residue modulo Min is a square. 

The following result due to Springer are interesting and important. 

Theorem 8 : Let R be a 2-Henselian discrete valuation ring in which 2 is a unit. 
Denote by k the residue class field R/M. Then the canonical suijection of R onto k 
mduces canonical isomorphism of their Witt rings, namely, s W{k). If K is the 
field of firactions of R, then W{K) = W(]c) © W{k) and there exists a ring isomorphism 
of W{K) with W(]c) {T\l{f-\). 
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Relationship between Quadratic Forms and Ordered Fields 

Let (K,P) be an ordered field and (Kq) be a quadratic space over K We say (l^q) 
IS positive definite if q(x) > 0 for all x e V, x 5 * 0, and negative definite if q(x) < 0 for 
all non-zero x eV, where the ordering is that induced by the positive cone P in K. 
Now given a quadratic space (f^q) over (K,P), we can decompose V as an orthogonal 
sum V"*' © V' where F^and V~ are such that V'^ with q restricted to it, is positive 
definite and F'with q restricted to it, is negative definite. (This orthogonal 
decomposition is independent of the choice of q from its equivalence class). The 
dimension of and F" are independent of the choice of orthogonal decompoistion 
of q and the dun. F^ -dim. F" is called the signature of the quadratic form q and this 
is an mvanant of q. 

The above result is the famous theorem of Jacobi and Sylvester and is known as 
Inertia theorem. Now using diagonahzation of quadratic forms, it is easy to*see that if 
< 3 ? and Y are two quadratic forms over (K,P), then sgn (ip J. ip) = Signature of the 
orthogonal sum of cp and ip = sgn (cp) + sgn (ip) and sgn (q) 0 \p) = sgn (cp) sgn(\p) 

Further, signature of the 1-dimensional form represented by q(a) = for all a e 
J^, n 0 is 1 and of the 0-form is 0. Thus this signature gives rise to a homomorphism 
of the Witt ring of K to the ring Z of integers. 

Suppose {K,F) is an Euclidean field so that every sum of squares is itself a square. 
Then fV(K) is isomorphic to Z and K* has exactly two square classes. In fact, for any 
formally real field AT, the additive group of its Witt ring fV(K) is not a torsion group 
(and this trivially implies that W(K) is not a 2-torsion group). 

As in the case of global fields where we have the local global principle of Hasse- 
Minkowski, for formally real fields we have Pfister's principle. We state this as the 
following • 

Theorem 9; Let AT be a formally real field. Then an element q in W(K) is a torsion 
element if and only if signature of q is zero for every ordering P e X(K)- the space of 
all orderings of K. 

If 9 is a torsion element of W{K) then its order is a power of 2 . 

The above result may be stated differently. Let P e XiK) and Kp denote the real 
closure of k with respect to P. The inclusion of K in Kp induces a homomorphism of 
fV(K) in fV(Kp) and Kp being real closed we have W(Kp) s Z under the signature map. 
Defining p from W(K) to itW{K^ where the product is taken overall P in X(k), we get 
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a (group) homomorphism and kernel of p is the 2-torsion part of fF(K) The kernel of 
p considered as a ring homomorphism is in fact the nil radical of fV(K) It is also the 
Jacobson radical of fy(K). Thus we have JTi(A!) = torsion subgroup of fF(K) is equal to 
the Jacobson radical of ff'iJO which is the intersection of all maximal ideals of fV(K) 
and it is also equal to the nil radical of W(K). (Nil radical is the intersection of all 
prime ideals of W(K) 

We have an immediate corollary from the above. If K is not formally real, then 
X(K) is empty so that the homomorphism p ; W(K) -^TtW{Kp) = 0 as p nins over an 
empty set Thus W(K) is a 2- primary group. 

We have a complete descnption of all orderings of an orderable field K via the 
prime ideals of the Witt ring of K. 

'^irst we start with a : 


Definition 4 ■ Let P be a prime ideal of the Witt ring W{K) of a field K. We define the 
characteristic of P as the characteristic of the quotient field of W(K)/p. 

The ideal I(K) of even dimensional quadratic forms over AT is the umque prime 
ideal of characteristic 2. Here is a complete description of all the prime ideals of 
W(K). 


Theorem 10 : (Leicht-Lorenz”, and Harrison). There is a one-to-one correspondence 
between the set XiK) of all orderings of K and the set of all prime ideals p of W{K) 
such that W{K)h^ is isomorphic to Z. 


One way is easy. Thus if P is an ordering of K, then the ideal arising out of the 
signature mapping Pp from W{K) -> Z is indeed a prime ideal. For the converse, let P 
be a prime ideal of W(K) such that the ideals 5'^^ ^ which is the kernel of the 


composition map W(K) -> W(Kp) -> 



where Kp is the real closure of K with 


respect to p. 


All these ideals are different — The prime ideals described by the above theorem 
are in fact the minimal prime ideals whereas the prime ideals of finite characteristic 
are the maximal prime ideals. 

We shall conclude our discussion with pre-orderings, variations that are 
compatible with pre-order and Witt rings. 
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Let r be a pre-order on the field K and Xt denote the set of all ordenngs P of K, 
which contain T. This set AV is a closed subset of the Boolean space X{K) and hence 
Xt is itself a Boolean space. We have a mapping from the Witt nng of A: to the ring 
C{X{K), Z) of mteger valued continuous functions from X{K) to Z and from the Witt 
ring Wj{K) of K relative to the pre-ordermg T. This ring Wt{K) is isomorphic to the 
quotient of W{K) by the ideal generated by the forms <1, -/) as t rmges oyer the non 
zero elements of T. Here by <1, - /> we mean the quadratic form (xf -tc|) In other 
words, we have JVj{K) = W{K)/ S W{K)-{\, -t), t e There is an imbedding theorem 
(due to Dubois) from W(K) to Z) where X(K) is the Boolean space with the 

Harrison topology) The kernel of the map is the nil radical For instance, if for K we 
take a formally real field and for T the pre-order given by sums of squares of elements 
(?i0) of K then the kernel of file homomorphism from W{K) to C{X(JC), Z) is exactly 
the ideal in W{K) generated by the forms (1, -a) as a range over sums of non zero 
squares of K. This kernel is in fact, the torsion subgroup of Wt{K). The quotient ring 
Wi{K) = WK/WfK for the case T = sum of squares of is called the reduced Witt ring 
K 


We have considered the compatibility of order and valuation on a field K. 
Likewise, we can talk of compatibility of pre-order and valuation We call a valuation 
V on a field iST to be compatible with a pre-ordering T if v is compatible with some 
ordering P of K with P 2 T. If this holds for all P in Xt, we say v is fully compatible 
with T. We set Rt to be the intersection of all valuation rings of K that are compatible 
with the pre-order T. Then Rt = R{P) where the intersection is taken over all 
ordenngs P in Xt For instance, if we take for T the sums of squares of non zero 
elements of K, then Rt is the real holomoiphy ring which is the intersection of the 
convex hull of the field of rational numbers with respect to all the orderings of K. 

We can associate another sub nng /?’ of K with each pre-order, namely the sub 
nng generated by the valuation rings R{P) associated with all the orderings P in Xp 
This R is indeed a valuation ring (being the over ring of valuation r ings ) and the 
valuation associated with this is the finest valuation of K which is fully compatible 
with the pre-ordering T. For more details and connections with ordered fields, etc., we 
refer the reader to the beautiful exposition of Lam* ’ . 

In conclusion I wish to thank Professor Ron Brown of the University of Hawai for 
his critical comments of an earlier version of this paper. Professor K.T. Arasu of 
Wright State University for useful discussion and Wright State University for their 
excellent facilities offered to me during my visit in 2001-02 and for the University 
Grants Commission in India for awarding me an Emeritus fellowship during the initial 
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stages of this work. I am grateful to the referee for pointing out several inaccuracies 
and typographic errors. 
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Abstract 

Two theorems are obtained to measure the degree of approximation from the Taylor transform 
and several other results including the result of Mohapatra and Chandra have been obtained. 

(Keywords . Taylor sums/holder metiic/degree of apjaroximation) 


Introduction 

Let C 27 t be the space of all 27t-periodic continuous functions / on [0, 2it] with 
Fourier series ; 


a 

/ ~ ^ + 2 (a„ cos nx + sin nx) , 

^ n^l 

Smgh'’^ defined the space by 

= ife C2n : l/(x) -/(y) l<Ka, (Ix-^'l)} 


( 1 ) 


( 2 ) 


and the norm || . || a,* by 

11/!U* = ll/]|c + sup A“*y(x,y) 


( 3 ) 
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where 


11/11,= sup l/'(x)l. 

(4) 

0^x<2n 



(5) 


and y) = 0. ©(f) and o*(0 are increasing fiinctions of t If o(lx-_>'|) < A |x-j^l“ and 
o*(|x-_y|) < , 0<p<a<l;^ and K being positive constants, then the space Ha 

IS give by : 

Ha = {/B C^^ : l/(x) -Ay ) ! < K lx-y|“, 0 < a < 1} (6) 


Let {iSn} be a sequence of partial sums of the given series ^ c„ , where S„ = cq + 


fp=Q 


Cl + ... + c„, then Taylor means of the sequence (S„} given by Hardy^, are 


T„=Za^S„irt = 0, 1 , 2 ,...). 


(7) 


Jt=0 


If (a^) is given by 


= <r<i, |re| < i), 


( 8 ) 


ihsn. the matrix is called the Taylor matrix. 


Miracle studied the Gibbs phenomena for Taylor means of Fourier series and 
Forbes studied the Lebesgue constants for regular Taylor sums. Prossdorf^ structured 
several interesting theorems on the degree of approximation defining the Holder 
metric. Chandra and Mohapatra^ studied some very interesting results, which include 
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the results of Prossdorf. Singh’’^ gave very general conditions to study the problems 
on approximation in generalized Holder metric. In this note, we take the definition of 
gener^ized Holder metric as given by Singh'’^ and try to generalize the following 
dieorem of Mohapatra and Chandra’ for Taylor operators. 

Theorem Let 0 < p < a < 1 , then for / e Ha, 

|| 7 ;(/)-/||, ( 9 ) 

where TJf.x) is the Taylor mean. 

MainResuits 

Theo) m 1: Let ©(t) defined in (2) be such that is increasing and ©(f)/?’ is 
decreasing. For 0 < r < Vi and 0 ^ p < t| < 1, we have 

{B(7c/n")}'^|log nr” (10) 
Theorem 2 : Let ©(t) defined in (2) be such that 

J m“’(d = 0 {H (/)} , (11) 

where H(t) > 0 be such that 

t 

J H(u)du = O^Hi})},t (12) 

0 

Then, for 0 < /• < Vi, 0 < p < r| 1 and / e Hw, we have 
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To prove above theorems, we shall require the following : 
Lemma 1 (Forbes^) ; For 0 < r < 1, lr0| < 1, 0 < t < n and /i, given by 


l-re" = he", 0 = tan‘ 


rsint 

J-rcos/ 


we have 


{(l-r)/h}” < exp (-knt^), k 


2(1 -rf 


{(l - r)/ h]" - exp 


2(1 -ry 




Lemma 2 (Miracle^) : Let r > 0 and 0 be as given in Lemma 1, then for some constant 
K, we have 


0 (0</<n/2). 

1-r 


Lemma 3 : Let 


<l>x (0 =/(x+0 +/(x-0 - 2f{x), 


where / e then 
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(19) 

and 

l4',(i')-'l>,(')|^‘to(i/|). (20) 

We shall use the following notations used by Mohapatra and Chandra’, as 

00 

L(«, r, t, 0) = '^a^s\nifi-\-\l'2)t (21) 

k=0 

and 

L(«, r,t,Q)= [{(l - r )/ sin {(« + 1 / 2) / + (w + 1)9}], (22) 


for 0 < r < 1, |r01 < 1 and (l-r e’^ = h e 0 < t < n. For completeness, we give the 
proof of (22). 

Using 

(i-i-rr 

- / ,W1 > 

fc=o (1 -/•<])) 


we have 


k=0 k=0 







k=0 


, 4. = e 


It 


it/2 


{i-re-T 


= Ime‘ 
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= Ime' 




1-re* 


he 


-tB 


- {{\-r)/hy^^ Im exp[{(«+l/2) t+{n + 1)6} i] 
= {{l-r)/hy^' sin {{n + 1/2) i + (« + 1)9>- 


Proof of Theorem 1\ 

Following Zygmund*, we have 

00 

T„ (f:x) = t.a^s,(x). 

fc =0 


We write 


i.(x) - T. if.x) -fix) - '> * 


where L(n, r, t, 6) is given in (22) 

1 '-w - '-w I 





= 

i + 


^0 

uln' ) 


= /l+/ 2 . 


(23) 


( 24 ) 


Since jl-r! ^ h, we have by (20) and (22) 



A NOTE ON THE DEGREE OF APPROXIMATION BY TAYLOR MEANS 


481 


nir. 


A i 




(25) 


By Lemma 3, (22) and (15), we have 

/j =0(l) I ^exp {-knt^)dt 


Hire 


= 0(l/«) f ^-^{exp (-Jtn^")}i/? 
1' ^ 


=0{«""-‘o)(7C/«'')}. 


(26) 


Next, we estimate I\ and h, using (19) instead of (20), we find that 


I\ - I\\ + I\2. 


(27) 


Since |l-r| < h and sin 0 < 9 for 0 < t < %ln and keepmg m view lemma 2 


/„=0 


©(jx- J |•|(/7 4•l/2)t + (n + l)^fe^ + 


rt 

~r) 


dt 


■o\ 


Ttin 

©(jx->'|)(9(«) J (l + ^?^ +r/(l-r))<il 


/„ = O M\x-y\)], 
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since ? for 0 < / < 7t/« and n is large. 

Further 

7,2=0 co|x- 7 |) I fylt)\L(n,r,t,Q\dt , 

'Kin 

Since |sin {(«+l/2) / + («+l)0)l 1 and jl-rj < h, owing to (22) 

n/n^ 1 

7,2=0 I - dt 

Kin ^ 

= o[co^x-y|)logn], (29) 

then 

7, = 0 [©(Ia: - log «] . (30) 

Since, {{\-r)/h}" < exp (-knt^) for 0 < t < 


72=0 (D(jx->'|)J -Lip,r,t,^)dt 

L Tt/w' ^ 


= 0 co(jx->^)— J r^^{-exp(-A7tt^)}fi(t 

Kin' 



(31) 

Now using 


=/r''C 

(32) 
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we have 


and 


Thus 


/j = O [{0(71/ {o(|x - >'|)log nf ] 


Sup{i-l.(..4=Suptt^^ 

Uc = ^)- /(^)l = - >1)1 


(33) 


(34) 


(35) 

(36) 


Combining above results, we get 

l7-;(/; x)- /(x|^ = o [{»{^ /»')r’ K")"' + '’"■'}]■ (37) 

This completes the proof of theorem 1. 

Proof of Theorem 2 follows analogously as the proof of Theorem 1 with slight 
changes, so we omit details. 

Corollary: 

Let r = 1/2 and tj = a in Theorem 1, and o(<) be as given in (2) such that alfylt is 
increasing and o(/)/7^ is decreasing, then we have 
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(38) 


for 0 < p < a < 1. 

Deduction 1 : If we put ©(0 = 0(r“), 0 < a < 1 and 0<p<a<l, inthe above 
corollary, then the theorem of Mohapatra and Chandra^ follows, 

a 1 

Deduction 2 : Again, putting p = 0, (o(t) = 0(t ),r= — , we get 


(39) 

Deduction 3 : If we put H{t) = 0 < a < 1 and 0 < p < a < 1, then 

!c(/;=t)-/(i| = o[«-''“-»'(lognf“ ] (40) 


and if H(t) = log (n/t), a = 1, 0 < r < 1/2 we get 

Deduction 4 ; If we put p = 0 in above deduction, we get 

and 


(41) 


(42) 







, a = 1, 0 < r < 1/3. 
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Abstract 

The effect of small perturbation in the Conolis and centrifugal forces on the location of 
libration point m a new kind of 3-body problem, has been studied when 0 < AT < 1 . 

(Keywords . 3-body problem/perturbation/Conolis force/centnfiigal force). 


Introduction 


A new kind of 3-body problem has been considered by Robe’, in which one body. 
Ml is a rigid sphencal shell filled with a homogeneous incompressible fluid of density 
pi The second one, Mz is a mass point outside the shell and M^ is a small solid sphere 
of density p 3 supposed to be moving inside the shell. The problem constitutes a model 
for the motion of an artificial earth satellite Ms inside another one Mi, Mz being the 
eartii Ms is assumed to be infinitesimal Mass of the three bodies are mi, mz and ms 
respectively. 

There are three forces acting on Ms- 

(i) The attraction of Mz 

(ii) The gravitational force exerted by the fluid of density p 


Fa = (4/3) nGpi ms Mi Ms 


where Mi Ms denote the distance of centres of the first and third body 
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(iii) The buoyancy force Fg = (4 / / Pj 


The total potential action at m3 is 

-Gm^ Ir^ +(4/3)7tGpi{l-(pi /pj) }r,^ 12 


= -Gm2lr2 +Kr^ /2 


where n = M, Mi, Mi, K = 4/3 Jtpifl-pi/ps) 

and G is assumed to be unity, 

Shrivastava and Garain^ used the small perturbation s and e' in Coriolis and 
centrifugal forces in the above problem. They considered M3 to be infinitesimal, the 
orbit to be circular and pi = p3 j.e. AT = 0 and obtained the equilibrium point [- ^ + 
E'p/(H-2p.), 0, 0], Shrivastava and Garain^ also studied the stability of libration point 
under the same assumption. They found that the range of stability increases for the 
point lying in one part ard decreases for the pomt lymg on the other part of the 
straight line-250 e + 1298' = 0. 

Here we shall apply the small perturbation in Coriolis and centrifugal forces in 
Robe's problran. We consider the orbit to be circular, pi ^ p3 and 0 < AT < 1 . 


The Location of Libration Point 
The equation of motion in Synodic system of co-ordinate are 
x~2y-x = dvldx 
y + 2x-y = dvfdy 


z~dvfdz 


(1) 
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where v = ]il[{x 2 -xf + ^ + - (Jell) [(jc-xi)^ + / + z^], 

|1 = W 2 /(Wi + m 2 ), Xi = -| 1 , Xi = 1-fi, 0 < }!< 1 . 

For unit of mass, we consider mi + m 2 = I and for the unit of distance, the distance 
betiveen centre of M\ and M 2 = I . 

vO < (4/3)?c Pi = (4/3) 73 Pi / {mi + mo) < (4/3) r. pi *._i 

= (4/3f n pi / (, 4 / 3 ) (r be the radius ofMA 

= !/r^ 

C < (4/3) % pi (i-pi/ps) < (i/r^) (l-pi/ps) 

0: - pi 

0 < /f < (p3 - pi) / (;“^p3) 

let (P3 - pi) / {r^pi) < 1 

i.e p 3 -pi</'V 3 

=> p3 < pi 

then 0 < ^ < 1 . 

Now we consider the perturbation in Coriolis and centrifiigal forces with the help of parameter 
a and p, the unperturbed value of both being unity. Thus we may write the eqn. (1) in the form 

X - lay - Px = dv/dx 

y + 2ax-^y = dvf^ 
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z = dvfdz (2) 

where a = 1 + 6, |Ej < < 1 
P = 1 + s’, |e’| < < 1 
This equation (2) is written as 

x~lay = dCljdx 
)' + lax = dOildy 

z = dClldz ( 3 ) 

where Q. = (p/2) (x^+/) -((k/l) [x - Xi)^ + / + z^]+\il{{{x ; -x)^ + / + 

From eqn (3), we get the Jacobi integral 

0 = x^ +_y^ +z^ -ICl + c = \\i{x,y,z,x,y,z) 

Libration points are obtained by assuming 

V* = ¥y= Vr 
i.e Q* = fly = Qj = 0 

fl* = Px + M.(l-p - x) / [(l-p - x)^ + - ^x+fi) = 0 (4) 

( 5 ) 


fly = Py - py / [(l-p - x)^ + y^ + - ky = 0 

= [li /{(l-li - xf + y^ + + i] z = 0 


( 6 ) 
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Fromeqn (6), we get tliat 

III / (I-^i - x)' + / + zY"" + k]:^0 
z = 0. 

Puttmg y = 0, z = 0 in eqn. (4), we have 

Px +M, / ( 1 -^ - xf - k(x + n) = 0 =/(x) (say) ( 7 ) 

AY = p(i-P) < 0 (vp>l) 
y(0) = ^[l/(l-H)'-^]>0 
V 0 < 1 - < 1 and 0 < < 1 

> 1 

=>l/(l-q,)^ > 1 

0 . 

Hence one root of the equation lies between -fi and 0. 

Now, we consider that x = - ^ + p 
Eqn (7) reduces to 

P(p-p) + |i/(l - p + p - p)^ - /t(p - fi - f^) = 0 

=> p^ (p-^:) + p2(-2p-pp + 2k) + p(2pP + P -^) + fi-fiP = 0 (8) 

Now, we use the small parameter method to find the value of p from eqn. (8) by 
putting p =0. 
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we have (p-A^) p(p^ - 2p + 1) = 0 
either p-^ = 0 or p = 0 or p^ -2p +1=0 
If p-jt = 0, then P = A:, which is impossible. 

If p^ - 2p +1 = 0, then p = 1, which is also impossible. 

So, p = 0(p) = Ci\i + + Csp^ +.... 

Putting this value of p in eqn (8), we have 

i^-k) +3cfc2\i* +3(ci^ +C3 +cjc|)p* +....} 

+ (- 2p - pp + 2^) p^ + 2 c]C 2P^ + (cjCs + Cj ) p'* 

+ (2pp+p-A:){(?ip+C2P^ +C3P^ + ....)+p-pp = 0 

Equating the co-efficient in the same power of p to be zero, we get 

1 - P + c, (P - it) = 0 

=i>ci = (p - 1)/ (p-it) = EV(p-it) 

Cii^-k) + 2 cip - c,^(2p+ 2k) 

=> C 2 = (-2e') / (p-^)^ (neglecting the higher power of s') 

2pC2 + C3 (p-it) = 0 
=> cj = 4e' / (p-it)^ 
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psV / (P-^) - (2e'p^) / + (4s'p^) / (P-/t)^ + 

= 8'p / (P-Ar) [ l-2p/(P-^) + / (p-A:)^ + ] 

= [E> / (P-it)] [ l+2p/(P--t)]-' 


= s'fi / (1 + s' + 2p, - A:) 

This value of p satisfies the eqn (8) 

:.x = -ji, + s'p / (1 + s' + 2p. - k). 

Let L denote the equilibrium point then 

L = i [ -p + 8'p/(l + e' + 2}x - k), 0, 0] 


When no perturbations act then E,e' = 0 and L = L (-p, 0, 0),which agrees with the 
Robe's result. 
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Abstract 

The Imear stability and first ordei normalization in planar restncted three body problem has 
been studied when more massive and mfinitesimal bodies are oblate spheroid under the 
condition that distances fiom center of mass of finite bodies to infimtesima! body are same 

(Keywords restncted problem of three bodics/massive body/mfinitesimal body/oblate 
spheroid/normalization/equatonal moment of mertia/polar moment of mertia/ 
Moser's theorem) 


Introduction 

Let us consider M], M 2 and M are three bodies with masses mi > m 2 » m ia 
which Ml and M are oblate spheroid. A and C denote the equatorial and polar 
moments of inertia of inertia of the body M respectively where as Ai and Ci for Mi. 
r,{i = 1, 2) be distanees between the centre of mass of bodies M, and centre of mass of 
body M. Let be the co-ordinates of the infuutesimal body at any instant under the 

barycentric system of co-ordinates The units are chosen as follows 

( 1 ) = 1, ji = m 2 /{mi + m 2 ) and l-p, = miKjni + m 2 ) 

(li) Gravitational constant = 1 . 

(iii) Distance between the centre of mass of primaries = 1 . 

Subbarao and Sharma’ considered the problem with one of the primaries as an oblate 
spheroid. It is of mterest to consider the oblateness of more massive primary and 
infinitesimal body 
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Equation of Motion 

Similar to Kondurar and Shinkarik' the equation of motion of our structure is 
wTitten as : 


£-<17 = 50 'ey 

i'-'-2£ = ^/£> (2) 

0= f/2 (/-:-/)•;- Ui 
0, = «j |l/r5)+ (^o)/|2r^}+ \A c? _ )f''l ■ 

U^=m,{yr,)+(Aa)/{2rff' 


<s={iC-A)/A} 


(-1 £ a < 1) 


< oi < 1) 

=(x-xy +y\r^ ={x-xj +/ ( 3 ) 

(*i'0) te,0) are co-ordinates of the primaries corresponding to barrycentre as 

origin. We have Xi = -p and Xj = 1-p 

0 = (l/2)(x^ +/)+0 -r)A +(li)/(/ 2 )+{(l-pMa}/n^ 

+ (lL4o)/(r3^)-^ {(l - p) (4a, )}/(2r,^) (4) 

And r^==(x+\if -i- (5) 

r^Hx~ l + ]if+y^ 


( 6 ) 
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Existence of Equilibrium Point 

For equilibrium position, dOJdx = 0, dQJdy = 0 

I e. X - {(l - ^i) (x + - {n(x - 1 + \i)}lrl - (3^a/2) [ {(l - ^i) 

(x + n)/r,'}-^(x-l + ^i)/r/}]-(34o,/2){(l-u)(x + ^i)}/r/ =0 (7) 

and T [ 1 - - v)h\ }“ W ''2 )" - p)/r,' ) }+ (^/r/ ) 

-^4a,(l-p)/(2n’)}=0 (8) 

For triangular libration point x 0 and _y 0 
from (8), we get. 

[1 - fl - n) }- (fi/ rl )- (3^a/2)fl - \i)l /-,* }+ (^i/ ri)- 

^yl,cri(l-M.)/(2r,’)}]=0 (9) 

Multiplying eqn. (9) by x and then subtracting from (7) and then putting, ri=r 2 = r, 
we have - ^Aall) - {A\a\/(\-[i)l2} = 0 (10) 

Now we calculate the nature of r with the help of Descart's rule of sign. For the 
real liberation point, r should be real and positive. 

The changes of sign of the equation depends on the sign of ct and Oi. 

Case I ; When ct and Oj both are positive. The sign of the above equation can be 
wntten as + . 

Here we get only one change of sign. Hence we get at most one positive real value 
of r . Again the equation is of odd degree. So only one +ve value of r exists. 

Case II : When o is +ve and crj is negative the sign of equation can be wntten 
as + - -+. 
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Here we get two changes of sign. 

.‘.We get at most two +ve real value of r. 

Case III ; When cr is -ve and Oi is +ve. The change of sign is written as + - + 

Here we get three changes of sign. We get at most three +ve real values of r 

Case IV ; When a and Oi both are negative. The change of sign is as follow 
+ - + +. 

We get two changes of sign. Hence we get at most two +ve real values of r. 

We frame a table given below : 


Case No 

Sign of cr 

Sign of Gi 

Number of 
changes of sign m 
the equation 

At most number 

of +ve value of r 

1 

+ 

+ 

1 

1 

2 

+ 


2 

2 

3 

- 


3 

3 

4 

- 

- 

2 

2 


Now we have 


X = (1/2) - n and / = - (1/4) [From (5) and ( 6 )] 

:.y = ±V(/-^-l/4) 

Hence triangular liberaticsi points are U = U [( 1 / 2 ) - pi, V(r^-l/ 4 )] 
By case 1, number ofZ,^ = 1; by case 2, number of Z, 4 = 2. 

By case 3, number of U = 3; by case 2, number of £ 4 = 2. 

And U = U [(1/2) -\i, -V(/^-l/4)] 

By case 1, number of Z-j = 1; by case 2, number of is = 2. 
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By case 3, number of I5 = 3; by case 2, number of £5 = 2. 

Puttmg r= \, libration points become the libration point of classical case. 

Linear Stability 

Since equilibriiun exists at £4 and £5 with co-ordinates 

[(1/2) - p, V(r^-l/4)] and [(1/2) - p, -V(r^-l/4)] respectively. 

Firstly, we calculate the stability for £4 with the help of Shrivastava and Garain^ 
as follows 

n^= 3/(4r*) [1 + 5Aal{2r^) + {^,a, (l-p)/(2r^)}] (11) 

a^= {+V-l/4)/(2r*)} [3(l-a)p+ {15/(2r*)} 

/l0(l-2p) -.4iai(l-p)] (12) 

(1/r') [(/•"-1/4) {3 + 15.4CT/(2r") + 15^i0,(l-p)/(2r")}] 

( )° indicates diat (x j') are replaced by (xo, yo) after proper differentiation when (xo ^yo) 
is the equilibrium point. 

.'. Characteristic equation regarding £4 can be written as 

X" + [4-(3/r^) {l+5^o/(2r*) + 5^,a,(l-p)/(2r^)} - 

{(r^-l/4)/r*} {3+15^a/(2r^) + ISyfiO, (l-p)/(2/'^)}] + 3/(4r’°) 

{1 + SAal{2r^) + 15>4iC, (l-p)/(2r^)} [(/'^-l/4)/r*} {3 + 15^CT/(2r^) 
+ 15^,01 (l-p)/(2r")}] - [{V(r^-l/4)/(2r*)} 3(l-2p) 

+{15^0 (l-2p)/(2r") - 15^,01 (l-p)/(2r^)] = 0 (14) 
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Eqn. (14) is a biquadratic equation in X. Hence it is a quadratic equation in 
{-Q±^{Q^-4PR)} / i2P) 

where, ? = 1 , 


Q = 4- {3/(5/) [{1 + 5Ja/(2/) + AiO^ (l-|a)/(2/)} {(/-1/4)//} 
{3 + 15^o/(2/) + IS^iC (l-n)/(2/)}] 
and R = 3/(4r’“) (1 + 5^a/(2/) + 5^ia, (l-ia)/(2/)} (3 + 15^a/(2/) 
+ 15^,01 (l-^)/(2/)} (/-1/4) - [V(/-l/4)/2/) {3(l-2n) 

+ 15^a/(l-2^)/(2/)- {15^,0, (1 -h)/(2/)}] 


For critical mass |ic, 2^ - 4PR = 0 and the characteristic roots are of conjugate pair. 
Hence I 4 is stable in this case otherwise L 4 is unstable. In the similar way the 
condition for stability of £ 5 ( 1 / 2 - 10 ., -V(/-l/4) can be obtained. 

First Order Normalization 

The Lagrangian function of the problem is written as 

L = ( 1 / 2 ) (jc^ + >'^ )+ (xy + yjc) + ( 1 / 2 ) + >^^ ) + (1 - [.t)/G + M /'a + 


pa)/ 2 } {(1 - \i)l + p/ }+ {(^a, )/ 2 } {(l - (r)/ 1 


where = (x + [0)^ + y'^,r^ =:{x-l - and triangular liberation point £4 is 
[(1/2) -p, ■'/(/-(1/4)}] in which r is obtained from the eqn. (10). 


£4 can be written as 


(r/2,r'/2). Now we shift the origin to £4. 


Thenx->x + (T/2),y-^>' 
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Eqn (15) reduces to 

L = (l/2)(x^ +/)+ + (Y/2)} j- {y + (Y'/2)} x] + (1/2) [{x+(Y/2)}^ 

+ {y + (Y'/2)}^] + (l-li)//-i + \i/r 2 + {(Ar^)l2} {(l-^i)/rj^ + |r/rf } 


+ {(^iOi/2} {(l~M')/^i } (^6) 

where 

r,-’= [{X + (1/2)}^ + {T+(Y'/2 )}^]-‘'^=/(x.y) (17) 

r;‘= [{x - (1/2)}^ + { 7 +(Y'/ 2 )}^ r'^^ = g(x,y) (18) 

[{X + (1/2)}' + {T+(Y'/2)}' ]''''=/;(^.T) (19) 

^ 2 '^= [{x - (1/2)}' + {>'+(Y'/2)}' Y^'^ = giix,y) (20) 


Expanding these in power series and substituting them in eqn. (16) and arranging the 
terms in ascending power of x, y, L can be written as £ = Lo + + 7-2 + 7-3 + 7,4 +.. and 

Lagranges equation of motions are 

i.e. x-ly= {l-(l/r')-3/la/(2/)-3^ia|(l-p)/(2r^)+3/(4r'^)+15 Acl(%r') 

+ ISyliO, (l-p)/(8r')} X + {3V(4/''-l)/(4r^) + 15^cTV(4r''-l)/(8r') 
- 3pV(4r'-l)/(2r^) -15.4cjpV(4r'-l)/(4r') 

+ 15.4|CT| (l-p)V(4r'-l)/(8/-')} y (21) 

and 

T + 2x = { (4r'- 1 )/(4r*)+3-( 1 //•')- 1 5 A V(4r'- 1 )/(8r')-3 p>/(4r'- 1 )/(2r') 
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-15Aop.V(4r^-l)/(4r’) +15AiC| (l-n)V(4r^-l)/(8r’)}x + 

- 3Acl{2r^) - S^ia, (l-^)/(2r^) + 3(4r^-l)/(4/-^) 

+ \5Aa (4r^-l)/(8r’)+15^ia| (1-^) (4r^-l)/(8r^)}j (22) 

The characteristic equation is given by 

X" + X' {4-(3/r') -15^a/(2r*) + 15 ^, 0 ,( 1 -^)/( 2 a-’)} 

+ {3/(4r'°)} {l+5.40/(2r^)+5.4,o, (l-n)/(2r")} [{3+ 15^ 0/(2/-') 

+ 15^,0, (l-n)/(2r')} (4r'-l)/4)}]-[{V(4r'-l)/(4r')} {3(l-2|i) 

+ 15(1-2^) Aa/{2r^) - 15^,0, (l-^)/(2r")}]^ = 0. (23) 


Let X = Wi, X 2 = iwt, X 3 = -iwi, X 4 = -Wa, 


Then 

W +^ 2 = {M3/r^)} - 15>t0/(2r*) - 15^|0| (l-n)/(2r*) 
wlw^= {3/(4/r'°)}{l+5^0/(2rV5^|O|(l-fi)/(2r")}[3+15A0/(2r^) 
+ 15^,Oi (l-n)/(2r^)} {(4r^-l)/4} - [{V(4r^-l)/(4r’)} {3(l-2ji) 
+ 15(l-2n)^0/(2r^) - 15.1i0, (l-^)/(2/-^)}]^ 

By applying Moser's theorem^ and following Garain and Shrivastawa^ we have 
Pi = l/2[l-(l+(5/3) (r*V) + (25/9) (r’^- 2r’-° + r’) - (64/225)r’°/(4/-^-l) 


(4-3/r^- 5r* + / [1+5 Aa/2r^) - (5/4) A^Oi/P] 

P 2 = l/2[l-{l+(5/3) (r*V) + (25/9) (r’^- 2r*“ + r") - (16/100)r"’ 
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2 

(4-3/r^ - 5r^ + 5r^) / (4r^-l )2 ]/[H-5 jG/2r^) - (5/4) J,a,/r^] 

The Hamiltonian corresponding to Lagrange's function is given by 

H=-L+ Py:X+ PyY. 

= 1/2(P/ +P;)+(yP. -xi>,)-(l-p)/n -p/r, -(^a/2){((l-p)/n^ +pr/) 
{(4a,/2}(l-p)/r,' 

Let us transform x x + yl2,y ^ x + •^12, P^-^ Px- yV2, Py-^Py + y/2, 

We get Hamiltonian 

Ho + Hi H2 + H3 + H4 + 

where, Ho = - (yVS) - (yV8) - (\/r) - (A(j/2r^) - AiOi (l-p)/(2r^) 

Hi=x {-(y/2) + l/(2r^) - p/r^ + 3^a/(4r^)-^^ 

2r 

- 3Jia, (l-p)/(4r^)+y'y{-l/2 + l/(2r^) + 3Aa/(4r^) 

-3.4,ai (1 -p)/(2a-')} 

H, =1/2{f^^ + ^/) + (yP^ - xPy)+x^ { 1/(2P) - 3/(8r^) + (3/4) Aa/r^ 

-(15/16) ^a/r’ + (3/4)^, a, (l-p)/(4r’) - (15/16) ^,a, (l-p)/r’} 
+y^{\l(2r^) - (3/8) (4r^ - l)/r^ + (3/4) ^a/r^ 

- (15/16) (4r^ - l)/r’+ (3/4) ^, 0 , (l-p)/(4r^) 



504 ASHOK KUMAR and D N. GARAIN 

- (15/16) ^,01 (4r'-l) (l-li)/r'} + y^xy{lr3IA)lr^ + (3/2)^/^' 

- (15/8) ^0//-’ + (15/4) ^0/r"- (15/8) ^,0, (l-^i)/r'} 

Now we find the canonical transformation from phase space (x, y, p^, py) in to 
phase space product of angle co-ordinates (<i)i, (t) 2 ) and action moment (/i, 1-^ and first 
order m , so that second order part of the Hamilotonian is normalized. For 

this we follow the method adopted by Whittaker^ 

We consider = dH^Jdx, -XPy = dH^dy, he = dH.JdP^, Xy = dH 2 /dPy 

which gives the following results . 

2bx + dy + XPx - apy = 0 

dx + Icy + apx + XPy = 0 

Xx - ay - Px = 0 

ax + Xy - Py = 0 (24) 


where, a = 1 

b - (1/2)//-^ - (3/8)/r^ + i3l4)Aalr^ - (15/16) Aa/r'^ 

+ (3/4) ^, 0 , (l-p)/r^ - (15/16) .1101 (l-M)//-’ 
c = (l/2)/r^ - (3/8) (4r^ -l)/r’ + (3l4)Ae5lr^ - (15/16) Aci{4r^-\)lr' 
+ (3/4) ^, 0 , (1-p)//-’ - (15/16) ^101 (l-p)(4r'-l)/r^ 

£/ = y' (-3/4)//-^ + (3/2)p/r* - (15/8) Aalr' + (15/4) Aa\ilr’ 

-(15/8) .4i0i (l-p,)/r^ 


solving eqn. (24), we get 

xl{-d + laX) = yl{2b - a^) = PJiaX^ - Xd - 2ab +a^) 
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= PyUX^ + a2X +2bX -ad) = K. 

we have foui values of A, and consequently four values of K. We consider the 
transformation (x, y, P^, P^ to (qi, q^, px, pi) similar to Garain and Shrivastava* and 
get 



where / = V(-i) 


Again we transform {qx, qz, px, pi) to (Qx, Q%, Px, Pi) similar to Garain and 
Shrivastava^ and get 

Qj = i2Ij/Wjy'^ sin «D„ Pj = i2IjWjy'^ cos (/ = 1, 2). 

Then the canonical transormation of phase space (X, Y, P^, P^) in to phase space 
product of angle co-ordinates (<I)i, O 2 ) and action momenta (Jx, li) is given be 


(X^ 


(QA 

Y 

— A 

Q2 

P. 

— /i 

P. 





where A = Jx J 2 =(Or*) (r, ^ = 1, 2, 3, 4) 

<3n = Xi - X3 iM'i/2 = kx(2aXx-d) -ksiWxOaXi - d)l2 

CI 12 = X2 - X4 ;m' 2/2 = k2i2aX2-d) -k4iW2i2aX4 - d)l2 
0 = 1, 2, 3, 4) can be made to satisfy the two condition Oxx = 012 = 0, which gives 
kx! ((m'i/ 2)(2 owi - id)} = k2l{2aiwx - d) = hx 
^ 2 ! {{'W2l2){2aw2 - id) = kj {-{2aiW2 - d)} = ^2 


other values of can be obtained 
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Ors - ^rs (k, Wl. 

where = 4 wJ +4-3/r^ -(l5/2)>la/r^ -(l5/2).4,a,(l-^)/r^ 

= 4W J + 12 - 3 / r ^ - 6Aa I -{9/ 2)Aiai (l - |i)/ r ^ 

=4w,'+3/r'-(9/2)^a/r'-(l5/2)^,ai(l-ii)//-', (j ^ I, 2) 

Thus H is normalized to 


H=WxIx-W 2 h 

Taking H = Ho + Hi + Ht, the equation of motion can be written as 
dl/dt = -dH / dQ>„ d<J>/dt = 0// / a/;, ( 7 = 1, 2) 
dl/dt = 0, d<t>i/dt = Wi, d<i> 2 ldt = -wj. 

The general solution is 

Ij = constant. (/ = 1, 2) 

O; = y^jt + constant. 
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Abstract 

The aim of this paper is to evaluate two integrals involving product of the //-function of two 
vanable" H [x, y] and the generalized polynomials S"' [x] • The values of the two 

mtegral jre obtained in terms of \j/ (z), the logarithmic denvative of f (z). 

(Keywords ■ //-function of two variables/generalized polynomials/logarithmic derivative). 

Introduction 

(a) The general class of polynomials 5*^’ "'^[x] is defined and represented* as 
follows: 


^“l. .»»r 

nj. .»r 


M = 


Inj/mj] 

I •• 

* 1=0 


["r/'Wrl r 

I n 


kr^O j=l 


~1J~ 




( 1 ) 


where n, = 0,1,2,...; rMi,...,»Jr are arbitrary positive integers, the coefficients 

A,,t, ^ 0) are arbitrary constants, real or complex. 


(b) The JT-function of two variables is defined and represented^ as follows: 



508 


PRAVEEN AGARWAL and C L. KOUL 


j\.P2 

= - <t>(^>n) 0i(^) e2(n) d^ (2) 


H[x,y] = H 


0,ni i«2,n2,m3.»3 


where 


f[ r(l-a^ + + 

(i,; = ^ (3) 

n TiOj - A^r]) f] r(l-6^ +p ^ + B,^) 

; = Ml+l ;=! ^ 

ft r(i-c^ + Y^4;f} 

0, iV = (1) 

^ JS- SL 

nns-Y;^) nrii-d^+bfi) 

J = »2+l j = nt2+\ 


and 02 (ri) is defined similarly. 

The contour L] is in the ^-plane and runs from — / oo to / oo, with loops, if 
neressary, to ensure that the poles off - 6, ^) (/ = 1, ..,/« 2 ) lie to the right and the 
po es of r (1- + Y, ^ ) (/ = l,...j^ 2 )j (i otj ^ T]) (/ = 1, ..,«i) to the left of 

the contour, Lz is similarly defined. 

It is assum^ that the existence and convergence conditions are satisfied for the 
various H~fancdons of 2-variables occurring here. 

By summing up the residues at the simple poles of the integral of (2), we have the 
following expansion^ for H [x, y]: 
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A=0 i=0 


hlH 


5} F, j 


( 5 ) 


where 


Cl =^,^ 2 =^.^.^= 0 . 1 . 2 , 

6 i Ai 


( 6 ) 


and 


/r^i,C2) = 4>(Ci,^2)e:(C,) e;(C2) 


( 7 ) 


where ^ (Ci.C 2) is given by ( 3 ), 


e; (CJ = 


ft r(l-c,+ y,CJ 



ft r(c,-y^) ft m-d^+SfiJ 

; = rt 2 +i ; = 


( 8 ) 


and ^2 (C2 ) ^ defined similarly. 


Ixt 11/ (z) 


x^iz) = 


r'jz) 

Tiz)- 


Main Integrals 

denote the logarithmic derivative of gamma function r(z) i.e. 


Also, let F(x)={2+Xi(1-x)+X2(1+x)} (9) 

and 


G(x) = 


(i-xr’a+x)*"' .-V 

{F(x)}""‘ *!■ 


^fr (i-yr’(i + ^)'' 
fi {F(X)}*'"'- 
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H 


z(i-xra+jc)''^ 

{£(x)}'‘2^''2 


( 10 ) 


We have 


J_j G(x)log 


(1-x) 

L£(x), 


1 [nrftny] 

dx = - y ... Y 


2 ^ 


- S SI 

^=0 kr^O A = 0 /=0 




-*-'i^-viCl-V2;2 


t=l 




Hog (1 + ^l)+W(a + s,k|+n^t;^+n^t;^)-^^f (a + b + k,(s,+t,) 


and 


+ (h +''l)Cl + (112 +''2)C2)}>'^' 


( 11 ) 


f G(x)log 


Hx) 

Eix) 


1 00 CO 

* = ^ L ■•• z z z 

*1-0 i.=Q li=o f=o 




M^-''1?1 -v 2?2 


8i^ 


{ Iogf7 + X2)+ ^fib + t,k, +v,^, +V2^2)~ 'f(a+b + k,(t, +5,) 


+ (h +V,)^j +(^j +V2)C2)}/^ 


( 12 ) 
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where 

= /(CnCz) Bia + s,k, +V.C, +V 2 C 2 X (13) 

/ = CbCs andy(?i,^ 2 ) are given by (6) and (7). 

(1 1) and (12) are valid for 

(i) \axgy\<^, \axg2\<Y> (14) 

where 

C' =- i - i Pj - I + 2 r, - S r, > 0, 

j^Ni+\ jr=l /=! ; = 1 y=i^2+l 

i - i fi. - £ F, + £ £, - i > 0 . 

;=J^i+l j = l j = l y=l y=V3+I 

(li) a, b, Hi, ji 2 , Vi, V 2 are all positive, and 


Re{a) + p, [i?e(c/i/5i)] + ii, [Re(f,/F,)] > 0 , 

Re(b) + v,[i?e(c/,/8,)] + v^ [Re(f,/F,)] > 0 . 

To establish the results in (1 1) and (12), we require the following integral : 



(l-x)‘’'^(l + x)*~^ y”*- 


cfr 

f=1 {Eix)r^‘’ 
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H 


O.A'i 

l\,Ql Pi.Q2<P3,Qi 


z(i-xr(i+^r 


dx 


1 [i^] [»r/m,] 

I jkl=0 Jt^=0 


{ n (1 + (1 + 


^ai\ri+2 I.isr2 .i.n-3 
+2,0+1 


>-(i+xir'*‘a+^2r''' ^ j\Px 

z(i+xi ri^2 (1+^2 r''2 e • (dy h.sa . (/^ > 1,23 


( 16 ) 


1 [«l^l ['ir/m/'J « « f r \ 

s - 1. Ez 

l * 1=0 kr=fi fe0^=0[i=l K|! 




h\i\ 




(1 + 31, ^ 

" 5iF, 



(17) 


with ^2 given by (6) and F (^i, ^ 2 ) given by (13) and 

P = il-a-s,k,; Pi.li2)>(l-*-^A;vi,V2),(Oj;a^,i4^)i,q 

Q = (!-«-*-*, (i-, +t ,) ; n, +v,,ii2 +V 2 ), 


where z = I,..,,?-. 

To evaluate 1, we substitute for '„"^[*] fro® £< 1 ^. (I), change the order of 

integration and summation, put the value of F [ x, y ] from (2), change the order of 
integration and integrate the x-integral with the help of Gradshteyn and Ryzhik^ 
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{b-a + X(t-a)+^(b-t)}“^P 

r(a)r(p) 
b-a r(a + p) 

a* b,^ (a) > 0, Re (P) > 0, 6 - o + X (/ - a) + n (i - /) 0 }, interpret the t|- 

integral in terms of the i7-function of two variables, we arrive at the right hand side of 
(16). 

Method of Proof : The result in (11) is established by taking the partial derivative 
of both sides of (17) with respect to a. (12) is similarly established by taking the 
partial denvative of (17) with respect to b. 

Special Cases 

In view of the presence of a large number of parameters in (1 1), (12), number of 
special cases can be obtained. However, we give here only one special case 

If, in (1 1), we take Xi = Xa = 0, it reduces to the following: 




7=1 


(l-x)^'(l + y)‘- 

(2)i,H 


H 


y(l-r)W(l + x)''‘ 
7(l-3C)^*Ml + y)^^ 

(2)W+V2 


log 


(1-x) 


dx 


2a^b-X 


I - S SI 

ki=Q kr=0 h=0 i=0 


n 

7 = 1 


k,\ h.UJ 5iFi 
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+ (K2+''2K2) 



( 19 ) 
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Abstract 

Making use of a known result we establish summation formulae for basic hypergeometric 
functions of two variables 

(Keywords : basic hypergeometnc fiinction/summation/Appel function) 


Introduction 

Snvastava and Jain' established the following identity. 


I 


/,m=0 


^i+mO^;q)i 


iq',q)i iq,q)m 


±Q„(Xp,q)„-^ ( 1 ) 

where is a bounded sequence of complex numbers and the parameters 

A, and p are essentially arbitrary 


Notations and Deflnitions 


A basic hypergeometric series is given by 



516 


PANKAJ SRIVASTAVA 




a\, at,. 

b\,bt , 


• . , Or', q ,z 
.■.,bs-, I 



n=Q 


(^1 ? ^2 ? ^ ^ ) n 

(q,bi,b2,bt, .. .b, , q)„ 


( 2 ) 


where 


f-1 

^2y 


_ n{n~l) 


and ipx, at, . , Or , q)„ ={a^,q)„{,at,q)„ {ar,q)„ (3) 


with the ^-shifted factorials defined by 


{a-,q)„ 


j 1 , if « =0 

[(l-ar)(l-o^) (l-a^"‘'),7/« = l, 2, .. 


{a;qU = 


I 


k=^0 


i\-aq^) 


( 5 ) 


(a ; q)>, = 


(g; g)« 

{aq”,q)a. 


( 6 ) 


For convergence of the senes (2), we need \q\ < \ and |z| < oo , when / = 1, 2, 
or max ( 1^|, |z|) < 1 when i = 0, provided that no zero appear in the denominator 

A generalized basic hypergeometnc senes of two vanables is defined as 


(o) ; (*) ; (6') ; q,x,y 

C ,D-,El[(c)-,{d),{d')-, i,j,k 


E 9 




X 
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A B b’ , 

Tt {ar,q) m+n n (br,q)„ n (br,q)„ x" y" 

!2 l ; £=2 (6) 

C D D , 

^ (f-r tq'im+n ^ i.^r ,q)m (d^ ,q)n i/i j <?)m (? j 9)n 

r=l r=l r=l 


The double senes converges absolutely for all bounded values of the complex 
arguments x and y when i.j, k e N and \q\<l and also when j =j = k= 0 provided that 
max ( 1^1 , Ix^ M) < 1 

Appell functions of first kmd is defined as ; 


b,b',c-x,y]= ^ X 

m=0 n=0 


m+rt ib,q)m jb' , q)„ x"” y” 
ic;q) m+n {q,q)m {q,q)n 


(7) 


Main Results 

Here we shall establish following suirmiation formulae ' 


(a) If we choose 


ja . q)n 
(aX HZ , q)„ 


We get 




[a, X, ^i,aXz, iiz, z] = 2^1 


a ,X^; z 
a X ^ z 


in (1), 


( 8 ) 


where <j)^[x, y] is the Appell function of first kmd defined m [Slater^ Chapter 
9(9 11)] 

Now summmg the 24>i senes on the right hand side of (8) by making use of the 
result [ Slater^ , Appendix IV, IV-2, p. 247] we get the summation formulae 


p; a Xz,^lz,z] = 


(X\iz ,az , q)„ 


(aX^z z; q)^ 


(9) 
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Taking %= l/|x m (9), we get 

(j)® \a ,\k ,az ,\Lz ,z] = \ 


( 10 ) 


Again takuig a = 1/X (i in (9), we get 


♦» [(^ rt- ; ^ , I* ; ^ ^ ^ . 1= 


( 11 ) 


(b) Choosing Q.„ = 


{a,q-[a ,-q-ia ,b,c ■,q)„ 

{4a , -4a , aq/b, aqlc, aqlX\i; q)„ 


z=aqlbcX [i 


in(l) 


and making use of the summation formulae [Slater^, Chapter III, eqn (3 3 1.3)], we 
get 


<!> 


^a,q 4a ,-q 4a , b ,c ;X,\i 
4a , -4a , aq/b, aq/c, aqlX\i ; - ; - 


aq/bcX ,aqlbcX\i 


_ {aq,aqlbc,aqlbXyL,aqlcX\x,-,q):„ 
{aqlb,aqlc,aq/X\i, aqlbcX\i\ q)„ 

Taking X = q''', ]i=q''' m (12), we get 


( 12 ) 


a, q4a,-q4a,b,c q~'' , q 
4a ,-4a ,aq/b,aqlc,aq^*'"^^ ; 


aq 


i+r 


aq 


l+r+i' 


be 


be 
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{aq,aqlbc, q) 

r+s 

{aq/b, aqlc, q) 

r+5 


As 6, c -» 00 in (12), we obtain 




a,q4a ,-q 4a ■, ; X, ; p, ; 

4a ,-4a ,aq/X\i;-;-; 2,2,1 


(aq, q)« 
(aq/X^ ; q)„ 


Again taking 6, c -> oo in (13), we get 




a, q 4a , -q 4a , q'‘ aq^*' , aq^^'*' 

4a , -4a , aq^*''*’ ; - ; - ; 


(13) 


(14) 


= (aq, q)r*s 

As X — >■ oo in (14), we get 


3<l>2 


a,q 4a ,-q 4a I 
4a , -4a ; 


-aq 

3 


(aq;q)<o 


Putting a = 1 in (16), we get the well known Euler identity 
I (-1)^ = iq ; qU 

r=Q 


(15) 


(16) 


( 17 ) 


Again taking a = -1 in (16), we get 
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s 


r=0 


1 + 

! + ?'• 


(~g i g )r r(3r-l)/2 

iq,‘l)r 


= i-q ; 


(18) 


(c) Choosing Q„ = m (1) and makmg use of the summation 

(aX \Lzq ; q)„ 

formulae [Agarwal^, Appendix II formulae 13 for c = 0 ], we get 
[a; X, p; u X p. zq, p 2 , z] = 


(Xpzg, azq , g)cc 
(aXpz^, z ; q)„ 


[l - (a + X p) z + a X p z] 


(19) 


Takmg X = 1/p in (19), we get 

; azq\ ^ ] = i 

(d) Choosing 


( 20 ) 


Cl - ,qq\)niq, q)„ 

i^q"' ; q)n izqt' ;qi) izq‘‘b^ ; qqx)„ (Xp ; q)„ 


in the identity (1) and making use of the summation formula [Agarwal^; Appendix II, 
formula 17) we get, 


^ (?»?'’ ig)i+m (q\ ;gi);+m {zq°^^ ',qq\)i+m 
/.»=o (Xp, z^"** ; q),^„ {zqt^^ ; q^),^„ 

(^ ,q)i i^ ;q)m (itz)'z"' 
i^q" qx ; qqx)un, (? ; q)i (q ; q)„ 


X 
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i\-z)Q.-zq‘‘ q\) 
Putting q\= q m (21), we get 


y (g. g*" ; q)l^n, (A. ; q)t (m. , g)m ^ 

1^=0 (X-n,, zq°^^ ; zq’^^)i^„ (zq^*^ ; 


(P^y 

(g ; g); (g ; g)« 


(a-zg°) (1-zg*) 
(1-z) {l-zq‘‘*'>) 


z = \ , (22) yields 


(g , g);H-m (>■; q)i (p , q)m (g°'^*''* ; q^)hm p' 


(X-li ; q)i^„ iq‘‘*‘’ ’ ; q^)u„ (q ; q)i (g ; q)„ q‘ 


= (l-g^-'Xl-g''^) ( 23 ) 

(2) Putting for ^ in the identity (1) and then taking zq for z and 

Making use of the result [Verma and Jain"*; (1.2)] we get the following summation 
formula : 


4,<» ; X, ^ ; X ^ ; M 9 ; -?]• I 
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^ (z , qU ^ (-Z -JiTii , q)^ ^ 

1 (-fi- P; ?)« (“V^’ ?)«> 


A number of other mterestmg results can also be scored 
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